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EXECUTIVE  SUMMARY 


This  report  presents  the  theoretical  basis  of  a  mathematical  model 
which  is  capable  of  characterizing  and  quantifying  the  hydrodynamic 
forces  and  yaw  moment  induced  by  a  current  on  a  moored  ship. 

In  Part  I,  the  behavior  of  the  flow  induced  by  a  current  on  a 
floating  vessel  (any  monohull  shape)  has  been  investigated  under  the 
assumptions  that  (i)  the  fluid  is  inviscid  and  incompressible,  (ii)  deep 
water  conditions  are  assumed,  and  (iii)  free  surface  effects  (such  as 
wave  making  and  scattering  of  possible  incident  waves)  are  neglected. 

All  possible  current  angles  relative  to  the  ship  centerplane  are  included. 

On  this  premise  the  necessary  mathematical  formulas  are  provided  herein 
for  predicting  the  longitudinal  and  transverse  forces  and  yaw  moment 
due  to  an  incident  current  for  any  given  hull  and  any  incident  angle. 

The  theoretical  model  is  constructed  with  application  and  further 
development  of  free-streamline  theory,  strip  theory  and  low-aspect-ratio 
wing  theory  for  the  overall  solution.  Various  physical  mechanisms 
have  been  elucidated  for  their  relative  contributions  to  the  resultant 
forces  and  yaw  moment  under  different  operating  conditions  including 
in  particular  the  different  ranges  of  incidence  angle  of  the  current. 

A  synopsis  of  this  idealized  theoretical  model  is  given  in  the 
Appendix  of  Part  I. 

In  Part  II  several  extensions  to  that  theory  are  presented  to  encompass 
more  general  conditions  that  may  prevail  in  practice,  such  as  sheared  flow, 
finite  water  depth,  and  wave-making  resistance.  These  extensions  include: 
the  first  order  correction  for  the  more  general  case  when  the  free 
stream  velocity  has  a  (weakly)  sheared  distribution  in  both  the  vertical 
and  horizontal  directions;  the  contribution  and  effect  of  the  rudder 
force;  wave  radiation;  and  the  effects  of  shallow  water. 

A  summary  of  the  complete  mathematical  model  is  presented  in  the 
Conclusions  and  Discussion  section  of  Part  II. 
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INTRODUCTION 


This  report  presents  new  research  results  achieved  under  the 
contract  N62474-82-C-8271  sponsored  by  the  Naval  Civil  Engineering 
Laboratory  (NCEL) .  The  primary  objective  of  the  study  is  to  assist 
NCEL  in  the  development  of  a  general  purpose  mooring  analysis 
capability.  This  capability  is  based  primarily  on  computer  models 
which  provide  solutions  for  the  static  and  dynamic  responses  of  an 
arbitrarily  moored  surface  vessel  subject  to  given  wind,  wave,  and 
current  loads. 

For  the  design  of  mooring  system  and  dynamic  positioning  (for 
localized  deep  sea  operation,  for  instance)  of  a  ship,  it  is  funda¬ 
mentally  important  to  attain  a  valid  estimate  of  the  external  forces 
and  moments  induced  on  the  vessel  by  the  motion  of  the  fluids  surrounding 
it.  There  are  three  kinds  of  fluid  mechanical  forces  and  moment  induced 
by  the  motion  of  the  surrounding  fluids — namely,  those  due  to  wind, 
wave  and  current  loadings.  The  methodology  acquired  by  Owens  and  Palo 
(1982)  at  NCEL  for  estimating  the  wind  load  has  proven  satisfactory. 

The  so-called  'drift  force'  induced  by  incident  waves  is  generally 
rather  insignificant  under  the  circumstances  of  NCEL's  specified  needs. 

In  sharp  contrast,  the  current  load  is  known  to  be  of  primary  importance, 
particularly  in  some  regions  where  surface  currents  up  to  3  to  5  knots 
(or  even  greater)  may  prevail.  But  the  existing  knowledge  about  current 
loads  and  the  very  few  experimental  results  available  (see  Remery  et  al., 
1973;  English  and  Wise,  1975,  and  as  reviewed  in  a  survey  article  by 
Palo  and  Owens,  1982)  are  insufficient  to  allow  them  to  be  reliably 
applied  to  wide  classes  of  ships.  We  agree  with  the  view  held  by  the 


NCEL  research  team  that  the  few  crude  methods  available  for  estimating 
the  current-induced  loads  are  inadequate  for  several  reasons.  First, 
for  those  methods  based  on  small-scale  experimental  results,  the 
information  is  insufficient  for  users  to  make  reliable  applications. 
Second,  these  methods  have  questionable  assumptions  that  appear  in 
conflict  with  physically  expected  scaling  behavior.  Further,  until 
recently  there  had  been  no  full-scale  tests  designed  for  validating 
these  methods.  A  series  of  full-scale  tests  was  carried  out  by  the 
NCEL  research  team  (as  has  been  discussed  by  Palo,  1982)  forming  the 
commencing  phase  of  an  important  task  (devoted  to  full-scale  tests) 
which  deserves  strong  support  for  further  development.  In  view  of  the 
urgent  need  of  conducting  theoretical  investigations  in  parallel  to  the 
NCEL  experimental  program,  the  present  project  has  been  undertaken  with 
efforts  concentrated  on  the  characterization  and  prediction  of  the  current 
loads  under  a  set  of  specified  operating  conditions. 

Current  force  is  a  hydrodynamic  force  which  depends  on  the  magnitude 
and  incidence  angle  of  the  water  stream  relative  to  the  moored  ship,  the 
characteristic  Reynolds  number  and  Froude  number,  as  well  as  the  geometric 
shape  of  the  ship  hull.  The  present  investigation  has  been  pursued  with 
the  primary  objective  of  providing  a  mathematical  model  and  methodology 
of  general  validity  and  practical  value. 

The  distribution  of  hydrodynamic  forces  induced  by  a  current  on  the 
vessel  can  be  divided  into  several  components: 

(i)  the  pressure  distribution  associated  with  the  wake-flow 
of  a  beam-on  current  leading  to  a  strong  transverse  form 
drag  which  in  turn  results  in  lift  and  drag  forces  and 


a  yaw  moment. 


(ii)  the  pressure  distribution  due  to  a  vortex-sheet  shedding 
flow  of  a  circulatory  nature  when  the  incidence  angle 
of  the  current  is  sufficiently  small, 

(iii)  the  frictional  force,  of  a  viscous  origin,  along  the  hull 
surface, 

(iv)  the  drag  contributed  by  the  various  appendages  such  as 
propellers,  rudder,  and  bulbous  bow  and  stern. 

Of  these  components  the  longitudinal  frictional  drag  is  generally  small 
compared  with  the  lateral  force  component  at  finite  incidence  angles. 

While  the  forces  contributed  by  various  appendages  can  be  considerable, 
they  are  not  emphasized  here  for  to  do  so  further  efforts  will  be  required. 

EQUATIONS  OF  MOTION 

A  ship,  which  is  free  to  move  as  a  rigid  body  except  for  its  being 

constrained  to  a  one-point  mooring,  has  six  degrees  of  freedom,  with 

three  in  translation  along  its  three  principal  axes  of  motion  and  three 

in  rota  :ion  about  these  three  axes.  It  is  convenient  to  employ  two 

cartesian  coordinate  systems  as  shown  in  Figure  1.  The  (x0»yo»z0)  system 

is  fixed  in  an  absolute  inertial  frame  of  reference  K  ,  with  the  origin  0 

o  o 

situated  at  the  point  of  mooring  (assumed  to  be  fixed  in  frame  K^),  and 
the  ZQ-axis  pointing  vertically  upward.  The  (x,y,z)  system  is  attached 
to  the  bod,  frame  of  reference  K,  with  the  origin  fixed  at  the  center  of 
gravity  of  the  ship.  The  current  velocity  in  frame  is  given,  denoted 
by  JJ  ,  and  is  assumed  (at  this  stage)  to  be  uniform  over  the  length  scale 
of  the  ship,  steadily  moving  parallel  to  the  horizontal  plane.  The 
prospects  thv  "he  jrrent  velocity  may  have  a  spatial  gradient  over 


4 


the  length  of  ship  dimension,  especially  when  it  forms  a  sheared  profile 
in  vertical  and  horizontal  directions  under  the  influence  of  shallow 
water  and  markedly  three-dimensional  bathymetry  can  be  admitted  by 
further  modification  of  the  method  presented  here.  Such  new  physical 
effects  can  be  significant  in  inducing  new  types  of  forces  and  moments; 
but  they  will  be  deferred  for  a  future  study.  Nevertheless,  temporal 
variations  of  the  current  velocity  can  be  admitted  on  the  assumption 
that  they  are  slowly  varying  when  compared  with  the  characteristic  time 
of  the  ship's  natural  oscillation  and  hence  can  be  regarded  as  quasi¬ 
steady. 

The  motion  of  the  ship  is  described  by  the  velocity  of  its  center 
of  gravity  and  the  angular  velocity  £1( t)  =  (0,0,d-y7dt)  about  the  vertical 
axis  through  the  center  of  gravity.  The  translational  velocity  of  the  ship 
relative  to  the  current  is  therefore 


(1) 


Its  x-  and  y-components  are 

u  =  u  -  u  ,  v  ■=  v  -  v  (2a) 

r  c  r  c 

respectively,  where 

u  =  U,  cos  0,  v  =  U,  sin  0,  u  =  U  cos  a,  v  =  U  sin  a,  (2b) 
b  b  c  c  c  c 


G  being  the  angle  of  inclination  of  the  velocity  vector  U  with  respect 

^  o 

to  the  x-axis  and  a  the  angle  of  incidence  of  the  current  velocity 
relative  to  the  ship's  centerplane  (see  Figure  1). 

The  equations  of  motion  for  the  moored  ship  are  conveniently  expressed 
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in  terms  of  the  body-fixed  system  of  coordinates  because  the  symmetry 
of  body  geometry  can  then  be  fully  exploited.  The  slow  oscillation  of 
a  moored  ship  can  be  described  as  consisting  of  three  modes  of  motion, 
namely  surge  (along  the  longitudinal  x-axis),  sway  (along  the  transverse 
y-axis)  and  yaw  (rotation  about  the  vertical  z-axis),  since  the  remaining 
three  modes  of  heaving,  pitching  and  rolling  are  usually  associated  with 
wave  excitation  and  are  relatively  insignificant  in  the  present  case. 
Accordingly,  the  equations  for  the  horizontal  motion  can  be  derived  from 
Kirchhoff's  equations  (see,  e.g..  Lamb,  1932,  Sect.  124)  together  with 
proper  modification  to  account  for  the  effects  of  the  prevailing  current. 
We  have  for 

surge:  (ws+  M, ,  )  -  (Ws  +  >n2a  )_Q  y  =  F*  +  Tx  ~  (  )  Cl  Vc  t  (3) 

sway:  (y)i5  +  yn22)^~-  -t  (  yns  +  mn  )£l  u.  —  Fy  +•  T^,  -  \xc  >  (4) 

yaw:  ^  +  (  ™2Z  "  )(U-UC  )  (V- Vc  )  =  Mz  +  ^  ,  (5) 

where  m  -  total  mass  of  the  ship, 
s 

mll’  m22  ~  X_  anc*  y^^P^ents  of  the  principal  added 
masses  of  the  ship, 

1  ,  I_„  -  moment  of  inertia  and  added  moment  of  inertia 

s  z  33 

(about  the  z-axis)  of  the  ship, 

F  ,  F  -  the  x-  and  y-hydrodynamic  forces  acting  on  the 
X  ■  ship  (in  addition  to  those  due  to  the  added  mass 
and  added  moment  of  inertia  of  the  ship) 

M  -  hydrodynamic  yaw  moment  about  the  z-axis, 
z 

T^,  -  x-  and  y-components  of  the  mooring  line  tension, 

-  distance  between  c.g.  and  the  mooring  hinge  of  the  ship. 


For  dealing  with  the  present  problem  of  slow  oscillations  of  ships 
moored  to  a  single  point,  in  the  absence  of  any  incident  waves,  the 
mass,  added  masses,  the  moment  of  inertia  and  the  added  moment  of 
inertia  can  all  be  assumed  to  be  constant  and  given.  We  proceed  next 
to  discuss  the  hydrodynamic  forces  and  their  origins  that  are  of  primary 
importance  to  the  dynamics  of  moored  ships. 


LONGITUDINAL  CURRENT  FORCE 

From  practical  experience  as  well  as  on  physical  grounds  it  is  known 
that  the  current  force  in  the  longitudinal  direction  is  generally  quite 
small  compared  with  that  in  the  transverse  direction  except  when  the 
current  runs  nearly  parallel  to  the  ship  length.  The  basic  reason  is 
because  the  typical  ship  form  is  slender.  Flow  separation  induced  by 
the  longitudinal  flow  would  occur  only  in  the  stern  region  whereas  viscous 
separation  of  the  transverse  flow  usually  occurs  along  the  entire  length 
of  the  leading  keel  line,  resulting  in  an  extensive  wake  formation  on 
the  downstream  lateral  side  and  hence  a  very  large  form-drag  coefficient. 

The  longitudinal  current  force  can  therefore  be  adequately  estimated 
by  the  conventional  viscous  resistance  formula: 

Fcx  =  i  L Ur  K)  cf  $ ,  (6) 

where  is  the  density  of  anbient  water, 

k  is  the  form  factor  (in  professional  use  for  given 
ship  form  class) , 

Cp  is  the  plank  viscous  resistance  coefficient. 


S  is  the  wetted  surface  area. 


For  all  practical  purposes  of  applications  to  hulls  under  commonly 
used  conditions  (the  surface  being  not  overly  eroded  or  rough)  we 
propose  to  adopt  the  1957  ITTC  line  of  C  (see  Figure  2),  namely 

r 

CF  =0.0TS(t>guRe-zj*  (7) 

where  Re  =  u^l/v  is  the  Reynolds  number  based  on  ship  length  l  and  the 
kinematic  viscosity  v  of  the  ambient  water.  This  viscous  resistance 
line  coincides  with  the  1956  ATTC  (American  Towing  Tank  Conference)  line 
for  Re  >  107  and  has  a  steeper  slope  than  the  latter  for  Re  <  107.  Thus 
our  choice  of  the  ITTC  line  should  be  conservative  for  prediction  of  the 
mooring  force  on  models  to  be  tested  and  for  extrapolation  to  the  prototype. 
While  there  is  no  universally  accepted  tables  for  the  form  factor  <  for 
various  types  of  ship  geometry  (including  the  generalized  kind  that 
allows  k  to  depend  on  both  the  Reynolds  and  Froude  numbers),  the  constant 
form  factors  k  which  have  been  adopted  in  current  professional  use  for 
nearly  all  categories  of  ship  forms  should  be  adequate  for  the  present 
purpose.  It  is  more  important  to  establish  a  method  capable  of  estimating 
the  transverse  current  load  with  comparable  accuracy  because  the  latter 
is  generally  much  greater  in  magnitude. 

TRANSVERSE  CURRENT  FORCE  AT  MODERATE  INCIDENCES 

The  flow  field  around  the  ship  generated  by  a  steady  current  at 
sufficiently  high  incidence  angles  (with  respect  to  the  ship's  center- 
plane)  resulting  in  a  viscous  wake  formation  is  very  complex  to  analyze. 
There  is  no  fully  adequate  theory  available  for  analytical  calculation 
or  numerical  computation  of  the  flow  for  Reynolds  numbers  lying  in  the 


range  encountered  in  general  practice.  For  a  best-estimate  engineering 
solution,  we  recommend  the  application  of  free-streamline  theory  to 
evaluate  the  wake  shape  and  form  drag  when  a  wake  is  formed  at  high 
enough  angles  of  incidence.  At  smaller  incidences,  the  oblique  wake 
on  the  transverse  side  will  cease  to  exist  and  we  adopt  slender-body 
theory  to  determine  the  hydrodynamic  forces.  It  is  therefore  convenient 
to  consider  the  two  regimes  of  (i)  small  incidence  angles,  without 
transverse  wake  formation,  and  (ii)  large  incidence  angles,  with  a 
transverse  wake  manifested. 

The  case  of  large  incidence  angles  can  be  most  effectively  treated 

by  means  of  strip  theory  of  free-streamline  flow.  The  wetted  hull  is 

divided  into  an  appropriate  number  of  strips  perpendicular  to  the  center- 

plane  of  the  vessel.  The  flow  in  each  strip  (except  for  the  end  sections, 

which  will  be  discussed  later)  is  taken  as  the  lower  half  of  the  two- 

dimensional  flow  of  a  double  body  which  is  formed  by  continuation  of 

the  flow  field  by  mirror  reflection  into  the  originally  undisturbed 

water  surface  (see  Figure  3).  The  flow,  after  encountering  the  leading 

side  of  the  hull  which  forms  a  blunt  surface,  is  assumed  to  detach  (the 

term  used  here  to  represent  approximately  the.  viscous  separation)  from 

the  body  surface.  The  detached  streamlines,  called  free-streamlines,  will 

envelop  a  near  wake  in  which  the  fluid  has  a  low  kinetic  energy  and  can 

be  approximated  by  a  region  of  constant  pressure  p  =  p  =  p100o  (at 

c  lou 

180°  point  of  the  base,  see  Figure  3).  The  force  so  determined  for  each 
strip  can  then  be  combined  to  yield  the  total  force  and  yaw  moment 
acting  on  the  vessel. 

To  solve  the  two-dimensional  free-streamline  problem  in  each  strip. 


we  adopt  the  general  theory  of  Wu  and  Wang  (1964)  which  provides  an 
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analytical  solution  in  terms  of  a  set  of  nonlinear  functional  equations 
for  computing  flow  variables  of  the  wake  flow  past  a  two-dimensional 
body  of  arbitrary  shape.  The  cross-sectional  shape  function  of  the 
'double  body'  may  be  prescribed  in  parametric  form  as 

y  =  $  (0<S<S)j  (8) 


where  s  is  the  arc  length  measured  along  the  wetted  portion  of  the  hull 
surface  (now  the  lower  half  of  the  double  body)  in  the  cross  section  of 
a  fixed  x.  Application  of  the  theory  of  Wu  and  Wang  to  the  present  case 
however,  will  still  require  some  extension.  The  detailed  procedure  is 
given  as  follows. 

The  set  of  nonlinear  functional  equations  can  be  expressed  as 


where  (see  Eq.  2) 
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In  these  equations  g(s)  denotes  the  angle  of  inclination  of  the  tangent 
to  the  solid  boundary  with  respect  to  the  coordinate  axes,  defined  by 


P(S)=  +»»  j  -ji-f-Ss-j  , 


here  expressed  in  terms  of  the  arc  length  s  along  the  boundary.  The 
variable  ;  is  a  parametric  complex  variable  in  which  plane  the  original 
flow  region  outside  the  near-wake  in  the  physical  plane  is  mapped 
conformally  into  a  lower  semi-circle  of  unit  radius,  the  point  5  (with 
|;j  <  1)  being  the  image  of  the  point  of  infinity  in  the  physical  plane. 

The  quantity  A  is  a  real  constant  depending  on  the  geometric  scale  of  the 
body  configuration.  The  transformation  (9)  and  (10)  maps  s(t)  onto  itself, 
i.e.,  s(0  is  a  fixed  point  of  the  transformation. 

This  nonlinear  functional  transformation  further  depends  on  the  end 
conditions  in  regard  to  the  point (s)  at  which  the  free-streamline  detaches 
from  the  solid  surface.  There  are  two  categories  of  detachment, 

(i)  fixed  detachment  and  (ii)  free  detachment  points. 

(i)  Fixed  detachment  —  The  detachment  points  of  this  category  are 
sharp  corners  on  the  body  surface,  from  which  free-streamlines  can 
separate  from  the  body  to  yield  a  valid  solution.  One  example  of  this 
case  is  the  cross  flow  past  a  ship  hull  equipped  with  bilge  keel  plates. 

In  this  case  both  detachment  points  (on  the  double  body)  are  known,  and 
equations  (9)  -  (13)  can  be  used  to  compute  the  solution  as  represented 


by  the  quantities: 


P  [(9(0;  <e;  . 
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An  iterative  algorithm  has  been  developed  for  numerical  computation  of 
the  above  parametric  solution.  In  nearly  all  the  cases  investigated, 
the  theoretical  results  given  by  this  model  have  been  very  satisfactory 
in  comparison  with  carefully  conducted  experiments,  such  as  the  lift  and 
drag  of  a  cavitating  hydrofoil  with  a  uniform  camber  of  16°  as  shown 
in  Figures  4  and  5. 

(ii)  Free  detachment  flows  —  By  free  detachment  we  mean  those  flews 
in  which  either  one  or  both  of  the  two  branches  of  free  streamlines  will 
separate  from  the  body  surface  at  a  point  where  the  surface  curvature  is 
finite  and  continuous  and  the  point (s)  of  detachment  is  unknown  a  priori 
and  must  be  determined  as  a  part  of  the  solution.  In  order  to  determine 
the  detachment  point (s),  we  must  impose  an  additional  condition  for  each 


free  detachment  point.  The  classical  condition  is 


A 
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as  -I, 

and/or  ^  — < ►  I, 


(15) 


as  the  case  may  be.  The  complexity  of  satisfying  this  end  condition  in 
the  numerical  solution  is  well  known  since  the  problem  is  highly  nonlinear 
and  its  solution  will  require  a  more  generalized  algorithm  and  a  new 
numerical  code.  Several  possible  approximate  methods  have  been  pursued 
under  the  present  program,  including  some  empirical  approaches  that  may 
effectively  aid  our  prediction  of  the  transverse  current  load  for  the 
case  when  flow  separates  from  a  smooth  keel  shoulder.  These  methods  will 
be  further  explored  when  the  development  of  a  numerical  code  is  carried  out 


STRIP  THEORY  AND  END  EFFECTS 


The  free-streamline  theory  described  above  will  enable  us  to  determine 


for  given  type  of  ship  form  (i.e.,  beam/draft  ratio,  block  coefficient, 
bilge  form,  etc.),  the  local  sectional  drag  coefficient,  Cpc(x;o)  say, 
due  to  a  beam-on  current  (i.e.,  at  sufficiently  large  angles  of 
incidence).  This  drag  coefficient  will  vary  with  the  longitudinal 
station  x  because  of  the  varying  sectional  hull  shape  and  will  depend 
on  a  nondimen sional  parameter:  i.e.,  the  wake-underpressure  coefficient, 

(fc>~hto’V  <’*> 

Here  is  the  pressure  at  upstream  infinity  and  P^qqo  the  wake  under¬ 
pressure  at  the  near-wake  base  (see  Figure  3) ,  both  with  the  hydrostatic 

part  deleted  (or  cancelled),  and  U  =  -v  =  v  -  v  is  the  speed  of  the 

yc  r  c  r 

beamwise  component  of  the  current  relative  to  the  local  ship  surface  (see 
Equation  2) .  In  engineering  applications  this  parameter  is  usually 
determined  by  experimental  means  due  to  lack  of  theoretical  recourse. 

Some  empirical  relations  will  be  sought  for  o  in  the  future  stage  when 
the  computer  code  will  be  developed. 

The  end  effect  on  wake  flow  past  the  ends  of  a  ship-like  slender 
body  is  an  unsolved  problem.  The  complexities  of  the  problem  lie  partly 
in  its  fully  three-dimensional  nature  and  partly  in  the  highly  variant 
curvature  of  the  hull  near  the  bow  and  stern.  It  makes  quite  a  marked 
difference  between  a  V-shaped  and  a  U-shaped  bow,  and  stern,  to  the  extent 
that  separates  the  fixed  and  free  detachment  categories  discussed  earlier. 
For  the  present  purpose  we  propose  to  apply  the  same  basic  free-streamline 
method  to  the  two  end  sections.  More  specifically,  each  end  section  will 
be  subdivided  into  an  optimum  number  of  horizontal  strips  together  with 
a  suitable  double-body  construction  in  order  to  apply  the  two-dimensional 


v.v.v.'jv’y *> vv.tj.v.v.v ’vwji'r; 7 ».« ■:« 77 *tt? ? ■••  tt?  .»  ■■■■  ■■?■?■ 
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free-streamline  theory  for  each  strip  of  the  end  section.  It  should 
be  qualified,  however,  that  the  validity  of  this  new  method  for 
correcting  the  end  effects  on  a  slender  wake  flow  will  require  further 
research,  with  both  theoretical  and  experimental  endeavors. 

The  total  transverse  current  force  can  then  be  obtained  by  direct 
integration  over  the  strips,  taking  into  account  the  variations  of  local 
current  relative  velocity  due  to  possible  yawing  of  the  ship.  Thus,  at 
moderate  angles  of  incidence  the  current  will  give  rise  to  a  transverse 
force  which  can  be  expressed  as 


r  i 


Fyc  “  [  <V*;°0(Uyc-nx) 


(17) 


-i< 


where  h(x)  is  the  sectional  draft  of  the  ship,  2.  and  £,  are  the  distances 

s  b 

from  the  wetted  stern  end  and  bow  end,  respectively,  to  the  center  of 

gravity  of  the  ship.  If  deemed  necessary,  the  magnitude  U  of  the 

yc 

transverse  current  velocity  may  be  allowed  to  vary  with  x.  Such  a  more 
general  case  can  arise  due  to  local  variations  of  bathymetry  under  some 
extreme  circumstances.  In  the  above  equation  (17),  the  transverse 
sectional  drag  coefficient  CDc(x;a)  also  symbolically  represents  the  end 
sections  for  which  alternative  options  with  further  subdivision  into 
horizontal  strips  will  be  explored  and  the  best  methodology  ascertained 
in  future  stages  of  this  program. 

Based  on  the  same  assumptions,  the  yawing  moment  acting  on  the  ship 
about  the  vertical  axis  fixed  at  the  center  of  gravity  induced  by  the 
oblique  current  in  this  range  of  incidence  is 

MC  =Z  fwj  CDC(^;^)(Uyc  -_Q*)  4 

-A 


(18) 


TRANSVERSE  CURRENT  FORCE  AT  SMALL  INCIDENCES 


When  the  ship  encounters  a  streaming  current  at  sufficiently  small 
incidence,  the  model  of  current-induced  wake  flow  and  the  strip  method 
presented  above  will  very  likely  cease  to  be  a  valid  model.  The  reason 
for  this  skepticism  is  the  argument  that,  on  physical  grounds,  such  a 
current  will  induce  vortex  sheets  to  be  shed  from  the  stern  (or  the  bow, 
whichever  is  downstream)  to  trail  behind  the  ship,  so  that  the  ship  will 
behave  more  like  the  lower  half  of  a  low-aspect-ratio  wing  rather  than  ?> 
stalled  slender  wing  as  delineated  in  the  previous  case  (with  a  transverse 
wake) .  We  therefore  assume  that  there  is  another  kind  of  hydrodynamic  force 
arising  from  a  circulatory  flow  origin  when  the  current  approaches  the 
ship  at  a  small  angle  of  incidence.  In  this  operating  condition,  we  shall 
regard  the  vertical  centerplane  of  the  double-body  as  a  lifting  surface 
of  low  aspect  ratio.  A  general  s lender-body  theory  has  been  developed  by 
following  the  approach  of  a  previous  theory  by  Wu  (1971)  for  analyzing 
arbitrary  motion  of  a  deformable  slender  body  (such  as  a  swimming  fish), 
equipped  with  side  fins  from  which  vortex  sheets  are  shed  as  a  result  of 
body  movement.  Now  with  the  contribution  from  an  oblique  current  included, 
the  theory  is  applicable  to  the  moored  ship  problem  at  hand,  especially 
when  longitudinal  vortices  are  shed  from  bilge  keel  plates  over  the  ship 
length.  As  these  vortex  sheets  are  convected  downstream  and  finally  enter 
the  stern  region,  they  can  have  a  significant  interaction  with  the  rudder, 
generating  lift  and  drag  forces  in  a  manner  quite  unlike  those  of  the 
wake-flow  origin  as  we  have  just  evaluated.  According  to  this  theory, 
the  lift  and  moment  on  the  ship  depends  on  whether  or  not  there  are 
vortex  sheets  being  shed  from  the  hull  surface  upstream  of  the  stern- 
rudder  section.  This  circulation  originated  flow  is  depicted  in  Fig.  6. 


.-a. 


In  the  case  when  there  are  vortex  sheets  being  shed  from  the 
bilge  keels,  we  shall  assume  that  the  shedding  starts  from  the  section 
at  x  =  physically  corresponding  to  the  section  of  maximum  transverse 
dimension,  and  continues  along  the  keel  plates  downstream  of  that  section 
where  the  keel  plates  become  trailing  edges.  In  the  presence  of  these 
trailing  vortex  sheets,  the  total  lift  in  the  direction  perpendicular  to 
the  incoming  current  (which  is  assumed  to  be  quasi-steady)  is 
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where 


Vxe  =-  uc  -  u  ,  Y CX,{)  =  Vyc  ~  n  *  ,  (20) 
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u  being  given  by  Eq.  2b,  U  by  Eq.  10b)  x  =  >,  at  the  stern-end 

yc  s 

.  is  the  station  where  vortex  sheets  ceases  to  Ko  shed  from  t:.e 
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and  „  A(x)  is  the  added  mass  of  the  transverse  segment  oi  the  chin  i 
w 

length  at  x,  for  relative  motions  in  the  >— uireet  io»j.  A(x)  •.•an  c 
calculated  for  given  ship  form.  We  note  that  the  added  mas:  ■  t  :e  .. 
is  half  the  added  mass  of  the  double  body  for  the  swav  n;  >de.  The 
drag  associated  with  the  lift  of  this  circulatory  f '  :>v  has  been  ..err 
from  the  theory  of  Wu  (1971)  by  assuming  quasi-stc-a.lv  iy;u  ■  xi:  .r  r  . 
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Transforming  the  force  svstem  (L  .  ,  I)  .  )  irto  t.h  •  •-•-  and,  \ 
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we  obtain 
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In  these  formulas,  the  stern  end  section  at  x  =  can  be  regarded  either 

as  the  ship's  stern  with  the  contributions  from  such  appendages  as  rudder 

and  propeller  units  superimposed  on  the  lift  and  drag  forces  or  as  an 

extension  to  the  trailing  edge  of  the  rudder  with  proper  account  taken 

of  the  effects  of  the  gap  between  stern  and  rudder.  The  details  of  these 

accounts,  however,  will  be  curtailed  here  but  can  be  accomplished  when 

the  development  of  a  computer  code  is  undertaken. 

A  direct  comparison  between  the  hydrodynamic  forces  and  yaw  moment 

in  the  two  operating  ranges  of  moderate  and  small  incidences  (see 

equations  17,  22,  and  23)  shows  that  when  SI 2,  <<  U  ,  the  transverse  force 

c 

F  of  the  wake-flow  origin  behaves  like  sin2 a  whereas  the  F  of 
yc  °  yc 

circulatory  origin  varies  linearly  in  sin  a  as  a-K).  In  order  to  achieve 
a  smooth  transition  between  the  two  theories  for  the  two  ranges  of 
current  incidences,  it  may  be  assumed  that  the  circulatory  flow  pre¬ 
dominates  in  a  range  of  small  I sinal  until  its  contribution  to  F 
becomes  equal  to  that  given  by  the  wake-flow  theory  and  drops  out  beyond 
that  transition  incidence  in  favor  of  the  latter.  This  assumption  is 


brought  forth  as  a  new  physical  concept,  which  seems  to  be  well  founded 
but  should  be  verified  by  experiments  for  its  full  validity. 


DISCUSSION 

The  theoretical  basis  for  a  mathematical  model  has  been  developed 
which  can  be  applied  to  calculate  the  hydrodynamic  forces  and  moment  on 
a  surface  ship  due  to  a  current  at  an  arbitrary  angle  with  respect  to 
the  ship.  Every  attempt  has  been  made  to  adopt  techniques  which  not  oniy 
represent  the  best  theoretical  approach  available,  but  also  can  be 
implemented  with  a  reasonable  computational  effort.  It  is  believed  that 
the  methodology  presented  herein  can  be  used  to  obtain  the  current  loads 
on  a  ship  moored  at  a  single  point  and  in  turn  to  deduce  its  resulting 
motion  subject  to  the  assumptions  of  (i)  deep  water,  (ii)  inviscid  and 
incompressible  fluid,  and  (iii)  no  free  surface  waves.  Even  though  much 
of  the  difficulties  due  to  the  highly  nonlinear  features  of  the  problem 
have  been  addressed  in  the  methodology,  there  are  still  a  few  remaining 
physical  factors  not  yet  included.  These  additional  issues,  which  are 
beyond  the  scope  of  the  present  project,  include  the  following: 

(i)  shallow  water  effects.  The  present  methodology  ignores  the  effects 

due  to  the  presence  of  the  ocean  bottom.  In  practice,  for  large 
ships  or  shallow  mooring  sites,  the  influence  of  the  bottom  can 
be  significant.  The  current  flow  under  these  circumstances 
would  tend  to  be  accelerated  in  passing  under  the  confined 
space  below  the  keel,  resulting  in  larger  transverse  loads. 

(ii)  spatial  variations  of  the  ambient  current.  The  ambient  current 

has  been  treated  herein  as  a  steady  uniform  one.  In  the 


ocean,  there  can  be  shear  currents,  i.e.,  current  velocity 


on  different  relative  importance  in  determining  the  current 
load.  Due  to  the  large  expense  of  field  tests,  they  will 
undoubtedly  remain  few  in  number.  However,  they  are 
indispensable.  Reliable  methods  for  extrapolating  model 
test  results  to  prototype  predictions  are  invaluable. 


RECOMMENDATIONS 


The  methodology  described  in  this  report  for  calculating  the  forces 
and  moment  due  to  a  current  on  a  ship  consists  of  several  separate 
formulations  each  applicable  for  portions  of  the  ship  hull  or  some  range 
of  the  parameters  such  as  angle  of  incidence  of  the  current.  The  overall 
application  involves  utilizing  all  these  formulations  and  integrating  the 
results  appropriately.  It  is  not  possible,  a  priori,  to  accurately  define 
the  exact  points  where  one  formulation  should  be  chosen  in  favor  of 
another.  This  must  await  the  actual  calculations  on  specific  hulls  and 
comparison  with  available  experimental  data.  It  is  therefore  recommended 
that  the  methods  described  herein  be  implemented  to  permit  the  numerical 
evaluation  of  the  forces  and  moment  induced  by  a  current  on  typical  ship 
hulls.  Computer  codes  should  first  be  developed  which  can  be  used  to 
calculate  the  forces  and  yaw  moment  acting  on  arbitrary  hull  shapes.  A 
set  of  representative  candidate  hulls  should  also  be  chosen  including 
ships  of  different  basic  designs  such  as  tankers,  destroyers,  cruisers 
and  so  on.  The  results  of  computations  for  these  hulls  under  a  variety 
of  conditions  should  then  be  critically  examined  not  only  for  purposes 
of  model  validation  and  comparison  with  experimental  data,  but  also  for 
insight  and  guidance  in  delineating  the  zones  of  applicability  of  the 
various  separate  formulations  so  that  a  unified  theory  can  be  achieved. 

In  addition  to  the  numerical  implementation  of  the  mathematical  model 
recommended  above,  effort  should  also  be  made  to  evaluate  the  effects  due 
to  the  various  physical  factors  discussed  previously  which  have  not  been 
incorporated  yet  in  the  present  formulations.  Some  of  these,  such  as  the 
shallow  water  effects,  can  be  of  great  practical  importance.  To  assure 


successful  applications  of  the  methodology,  it  is  also  essential  to 
secure  valid  scaling  relations  between  small  scale  model  test  results 
and  the  prototype  behavior. 
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NOMENCLATURE 


A  A  real  constant  representing  the  geometric  scale  of 

the  wetted  body  surface,  used  in  Eq.  9 

A(x)  A(x)  times  the  water  density  o  gives  the  added  mass 

of  the  transverse  segment  of  tKe  ship  hull  per  unit 
length  at  x 

b  Beam  of  the  ship 


C_  Plank  viscous  resistance  coefficient 

F 

C^^  The  transverse  sectional  drag  coefficient 

D  ^  Induced  drag  associated  with  the  circulatory  flow 

for  small  incident  current  angle 


F  Longitudinal  current  force 

cx  6 


F 

x 


The  x-  and  y-  hydrodynamic  forces  acting  on  the  ship 
(in  addition  to  those  due  to  the  added  masses  and 
added  moment  of  inertia  of  the  ship 
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The  x-  and  y-  hydrodynamic  forces  on  the  ship  due 
to  the  drag  and  lift  of  the  circulatory  flow 
induced  by  the  current  at  small  incidence  angles 
to  the  ship 

Sectional  draft  of  the  ship 

Moment  of  inertia  and  added  moment  of  inertia  (about 
the  z  axis)  of  the  ship 

The  form  factor  of  a  given  ship  hull 

Total  lift  in  the  direction  perpendicular  to  the  mean 
incoming  current  due  to  the  sheeding  of  vortex 
sheets  (see  Eq.  19) 

Ship  length 

Distance  between  the  mooring  hinge  and  the  center 
of  gravity  of  the  ship 

Distance  in  x  direction  representing  the  locations  of 
bow  and  stern  ends  to  the  center  of  gravity  of  the  ship 

Location  (in  x-direction)  of  the  maximum  transverse 
dimension  of  the  ship 

Location  in  x-direction  where  vortex  sheets  cease  to 
be  shed  from  the  hull 

Total  mass  of  the  ship 

x-  and  y-  components  of  the  principal  added  masses 
of  the  ship 

Hydrodynamic  yaw  moment  about  the  z-axis 

Yaw  moment  acting  on  the  ship  about  the  vertical 
z-axis  (see  Eq.  18)  due  to  the  wake  forming 
oblique  current 

Moment  of  forces  about  the  vertical  z-axis  due  to 
the  sheeding  of  vortex  sheets 

Pressure  at  upstream  infinity 

Wake  underpressure  at  the  near-wake  base 

Reynolds  number  defined  as  Ur*£/v 

Reynolds  number  defined  as  U^'b/v 


S 


Wetted  surface  area  of  the  ship 


s 


Arc  length  used  in  prescribing  the  cross  sectional 
shape  function  of  the  "double  body" 

x-  and  y-  components  of  the  mooring  line  tension 


The  x-  and  y-  velocity  components  of  the  ship's 
center  of  gravity 

U  =  (U  cos  a,  U  sin  a)  =  (u  ,v  )  the  x-  and  y-  components  of  the 

current  velocity 

The  x-  and  y-  velocity  components  of  the  ship  relative 
to  the  current 

v  The  beamwise  velocity  component  of  the  current  relativ 
to  th  e  local  ship  surface 

U  =  u  -u  The  longitudinal  velocity  component  of  the  current 

relative  to  the  local  ship  surface 

a  Incidence  angle  of  the  current  with  respect  to  the 

ship  centerplane 

Angle  of  inclination  of  the  boundary  slope  of  the 
hull  surface  in  x = constant  plane 

Two  nonlinear  functionals  defined  in  Eqs.  (11)  and  (12) 


A  complex  variable  forming  a  parametric  plane 
v*  A  function  defined  in  Eq .  (l.lb) 

v  Kinematic  viscosity  of  the  ambient  water 

p  Density  of  the  ambient  water 

w 

o  The  wake-underpressure  coefficient,  or  the 

cavitation  number 

The  angular  velocity  of  the  ship  about  the  z-axis 
through  the  ship's  center  of  gravity 


6(s) 

£[e] ) 


V.(t) 
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APPENDIX 


SYNOPSIS  OF  THE  PROPOSED  THEORETICAL  MODEL 


In  order  to  assist  the  principal  engineers  concerned  with  the 
continued  planning  and  accomplishment  of  the  project,  and  to  aid  the 
future  users  of  the  theoretical  model  which  is  being  developed  here, 
we  present  in  this  appendix  a  brief  orderly  outline  affording  a  quick 
general  view  of  the  physical  background,  simplifying  assumptions, 
conditional  qualifications  and  estimated  range  of  validity  of  the 
proposed  theoretical  model. 


These  equations  are  expressed  with  respect  to  the  body  frame  of  reference 
(fixed  with  the  moored  ship) ;  they  are  derived  especially  for  the  present 
case  in  which  we  have  a  current  of  arbitrary  magnitude  and  incidence 
angle  passing  about  the  moored  ship. 

Under  the  same  premise,  the  added  mass  coefficients  m^ ,  m^  anc^  ^33 
can,  quite  accurately,  assume  their  values  calculated  from  the  "double-body 
model,"  meaning  the  extended  flow  obtained  by  reflection  of  the  original 
into  the  initially  undisturbed  water  surface.  By  this  approximation,  m  , 
m^2  arid  I  attain  values  which  are  one-half  of  their  corresponding 
coefficients  of  the  double-body  flow,  which  is  unbounded  in  extent  and 
symmetric  about  the  mean  water  surface.  Clearly,  the  values  of  these 
coefficients  depend  only  on  the  hull  shape,  but  not  on  the  time  history 
or  the  frequency  of  body  motion. 

2 .  Longitudinal  current  force 

For  the  determination  of  the  longitudinal  viscous  resistance  induced 
by  a  prevailing  current,  the  formula  (6)  is  proposed  for  the  following 
reasons.  First,  when  the  current  runs  nearly  parallel  to  the  ship  length, 
the  problem  becomes  that  of  a  ship  moving  under  tow  or  by  self-propulsion, 
in  which  case  the  formula  is  well  based  on  the  abundant  towing-tank  test 
results  and  prototype  experiments  acquired  in  the  profession  of  naval 
architecture.  While  it  is  true  that  in  the  case  of  an  oblique  current 
past  a  ship  no  assessment  of  the  formula  (6)  has  been  made  available  (to 
our  knowledge),  the  existing  experimental  data,  though  very  scanty,  have 
confirmed  our  expectation  that  In  such  cases  the  longitudinal  viscous 
resistance  is  quite  insignificant  in  comparison  with  the  transverse  current 
force  when  the  flow  has  a  considerable  length  of  separation  along  the  hull. 


26 


3 .  Transverse  current  force  at  moderate  incidences 

Ln  this  case  large  values  of  transverse  current  force  coefficient 
(of  order  unity)  are  known  to  exist.  To  determine  this  coefficient  by 
numerical  methods  using,  a  fully  viscous  fluid  model  (e.g.  the  Navier-Stokes 
equations  or  a  turbulent  boundary  layer  model)  is  still  beyond  the  present 
state  of  the  art.  For  engineering  applications,  such  as  for  our  problem 
at  hand,  the  free-streamline  theory  proposed  here  is  considered  to  be  th^ 
most  effective  to  provide  a  two-parameter  family  of  solutions.  This  method 
is  different  from  the  singularity  method  that  employs  mass  sources, 
doublets  and  vortices,  the  expansion  method  that  uses  orthogonal  functions, 
or  the  boundary  integral  equation  method  based  or.  Freon's  integral  formula. 
The  great  power  of  this  method  is  iully  exploited  in  the  case  of  two- 
dimensional  plane  flow  by  means  of  generalized  conformal  transformation 
(involving  mathematical  branch  points)  and  complex  function  theory. 
Application  of  the  method  to  the  present  problem  of  three-dimensional 
ship  motion  is  effected  by  the  concept  of  strip  theory  which  appears  to 
be  applicable  to  hull  shapes  of  conventional  form.  The  idea  of  introducing 
the  strip  approximation  for  the  bow  and  stern  end  sections  is  rather 
innovative  and  should  require  further  studies  for  its  full  verification. 

It  should  again  be  noted  that  the  two  parameters  the  solvit  ion  will 
depend  upon  are  (i)  the  wake  under-pressure  coefficient,  and  (id)  the 
location  of  the  detachment  point  from  which  the  free  streamline  will 
separate  from  the  hull  surface  to  generate  i  cross-1' i :  wake  formation. 
Although  a  closure  condition  has  been  proposed,  in  terms  of  equation  (IV) 
of  the  text,  for  the  determination  of  parameter  (ii),  various  alternative 
closure  conditions  for  the  two  parameters  have  been  explored  in  the  present 


investigation,  including  several  empirical  approaches,  with  intent  to 
improve  the  accuracy  of  the  solution.  Furthermore,  extensive  efforts 
have  been  spent  during  this  study  in  examining  the  prospects  of 
convergence  of  the  numerical  iteration  of  the  nonlinear  integral 
equations  (eqs.  9-12  in  the  text).  This  concentrated  effort  was 
deemed  necessary  for  the  present  application  (to  hull  forms  whose  cross- 
sectional  shape  may  be  regarded  as  blunt  and  with  highly  curved  corners) 
because  our  experience  of  the  previous  success  of  the  theory  was  limited 
to  thin  bodies  with  sharp  corners. 

4.  Transverse  current  force  at  small  incidences 

The  slender-lifting-body  theory  proposed  here  is  believed  to  be 
particularly  well  suited  for  the  case  when  the  incidence  held  by  the 
approaching  current  is  small.  In  this  regime  the  cross  flow  is  expected 
to  become  so  weak  that  the  cross-flow  wake  formation  should  be  evanescing 
but  with  vortex  sheets  continually  being  shed,  only  now  confined  more  to 
the  stern  region.  Since  the  original  theory  cited  in  the  text  is  best 
applicable  to  thin  bodies,  such  as  thin  ships  of  the  Wigley  form,  further 
modifications  of  the  theory  have  been  explored  in  order  to  extend  its 
validity  to  cover  hull  forms  of  conventional  categories.  The  optimum 
modification  remains  to  be  identified  by  numerical  calculations  to  be 


carried  out  later. 
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Figure  4.  Comparison  between  theory  (Wu  and  Wang,  1964)  and  experiment 

(Parkin,  1958,  J.  Ship  Research  1.,  34-56)  for  the  lift  coefficient  C 
of  a  cavitating  circular  arc  hydrofoil  (of  semiarc  angle  y  =  8°)  and 
its  variation  with  cavitation  number  (or  wake  underpressure 
coefficient)  a. 


1.  Introduction 


In  Report  AMT  R-101  entitled  "Investigation  of  Current  Loads  on 
Moored  Vessels",  a  methodology  was  presented  with  a  proposed  theo¬ 
retical  model  for  calculation  of  hydrodynamic  forces  on  a  moored  vessel 
in  the  presence  of  a  uniform  current  in  deep  water  without  considering 
the  effect  of  surface  waves.  The  theoretical  model  was  constructed  with 
application  and  further  development  of  free- streamline  theory,  strip 
theory  and  low-aspect-ratio  wing  theory  for  the  overall  solution. 

In  this  present  sequel  to  the  first  original  report  (AMT  R-101), 
several  extensions  to  that  theory  are  presented  to  encompass  more  general 
conditions  that  may  prevail  in  practice,  such  as  sheared  flow,  finite  water 
depth,  and  wave-making  resistance. 

In  Sections  2  and  3,  a  methodology  is  formulated  in  detail  to  obtain 
the  first  order  correction  to  the  original  freestreamline  -  strip  theory 
model  for  the  more  general  case  when  the  free  stream  velocity  has  a 
sheared  distribution  in  both  the  vertical  and  horizontal  directions.  It  is 
based  on  a  small  perturbation  approach  and  is  applicable  to  the  case  of 
weak  shear.  Section  4  presents  results  for  hydrodynamic  forces  and 
moments  on  the  moored  ship  as  a  result  of  a  current.  The  particular 
contribution  and  effect  of  the  rudder  force  is  evaluated  in  Section  5.  Wave 
radiation  from  a  ship  moored  in  a  steady  current  while  not  a  major 
contributor  to  the  hydrodynamic  forces,  is  analyzed  in  Section  6.  Finally, 
Section  7  discusses  several  aspects  of  the  effects  of  shallow  water. 
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where  U  is  a  maximum  measure  of  the  current  speed  and  1  is  the 

c  max 

ship  length.  Such  typical  variations  can  generally  be  ascribed  to  local 
topographical  distribution  of  water  depth  and  changes  in  forcing  distur¬ 
bances  such  as  tides  and  distant  storms.  Variations  of  this  weak  category 
can  usually  be  approximated  by  its  Taylor  series  about  the  mean  water 
suface  (at  z  =  0), 

9u  9v 

uc  =*c(x.o,t)  +  ,  vc  =  vc(x,  0,  t)  +  z(_ -p)  .  (2) 

o  o 

Vertical  variations  of  the  current  velocity  beyond  the  terms  linear  in  z 
are  less  important  and  can  be  neglected  when  the  local  water  depth  hQ 
is  large  compared  to  the  ship  draft  d  (say  h^  >  2d).  In  very  shallow 
water,  a  steadily  running  current  has  its  velocity  varying  with  depth 
according  to  a  one-seventh  power  distribution  typical  of  turbulent  flow 
over  a  smooth  bottom.  If  the  clearance  beneath  the  ship  is  of  order  of  the 
ship's  draft  or  shallower,  the  effective  current  velocity  experienced  by 
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the  ship  may  require  terms  nonlinear  in  z  to  be  also  included  in  the 
expression  in  (2),  but  such  exceptional  cases  must  be  treated  separately 
since  the  shallow- water  effects  then  becomes  of  primary  importance, 
dominating  over  the  effects  due  to  nonuniformities  in  the  current.  This 
extreme  case  will  be  analyzed  in  a  later  section  (see  Section  7). 

The  first  terms  of  the  expansions  in  (2)  represent  the  velocity 
components  of  the  surface  current,  which  may  vary  with  time,  but  only 
very  slowly,  and  may  vary  with  x  along  the  ship,  but  also  very  gradually. 
Further  expansion  of  these  surface  current  velocity  is  unnecessary  since 
in  strip  theory  for  slender  ship  motions  the  weak  variation  of  current 
velocity  with  x  and  t  can  be  taken  as  being  parametrical  in  nature. 
Consequently,  the  strip  theory  developed  earlier  in  this  report  can  be 
directly  applied  without  modification.  We  now  proceed  to  consider  the 
problem  of  the  transverse  wake  flow  in  a  vertically  sheared  stream  for 
each  strip  of  the  ship. 


3.  Two-dimensional  wake  flow  in  a  vertically  sheared  stream 

We  consider  the  wake  flow  of  a  uniformly  sheared  free  stream 
past  a  two-dimensional  blunt  body  of  arbitrary  form,  with  the  stagnation 
streamlines  separating  from  body  surface  to  become  free  streamlines 
enveloping  a  finite  cavity  (or  wake)  of  constant  pressure,  pc>  Viscosity 
is  neglected  so  that  the  effects  of  the  vorticity  field  alone  on  free- 


streamline  wakes  are  considered.  The  basic  flow  has  the  velocity 


components  (0,U(1  f  yz),  0)  in  the  Cartesian  coordinates  shown  in  Fig.  1 
where  U  and  y  are  constant,  and  hence  the  flow  has  a  uniform  vorticity 
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It  follows  from  Helmholtz's  vorticity  theorem  that  the  resulting  velocity 
field  u(x)  of  the  wake  flow  will  be  characterized  by  the  same  uniform 
vorticity  distribution  throughout, 

<t  =  Vxu  =  =  -  Uv  e  x  ,  (4) 

where  je  is  the  unit  vector  pointing  out  of  the  paper  and  w  is  the 
vorticity  vector.  This  is  an  immediate  consequence  of  the  Helmholtz 
wrticity  equation, 

dco 

-at  =  {ar+  *  ‘  V)*  =  0  '  (5) 

which  holds  for  two-dimensional  motion  of  an  incompressible,  inviscid 
fluid  as  is  the  case  at  hand. 

Accordingly,  the  resulting  velocity  field  u  =  (0,  v,  w)  will  consist 
of  two  parts,  one  retaining  the  original  uniform  vorticity  throughout  and 
the  other  being  ir rotational,  hence  possessing  a  scalar  velocity  potential 
<j)( y,  z;y)  which  depends  on  the  vorticity  strength  -y  as  a  parameter.  Thus 
the  velocity  components  can  be  expressed  as 


TT  ^90  9* 

u  =u^z  +_Z  =  — 


90 


(8) 


In  addition  to  these  equations  for  the  kinematic  variables,  we  have,  for  the 
dynamics,  the  Bernoulli  equation 


P  + 


2  P<U 


2  +  w2)  =  H(tf) 


(9) 


where  p(x)  is  the  pressure,  p  the  (constant)  fluid  density,  and  H('I') 
represents  a  constant  along  each  stream  line.  In  terms  of  the  kinematic 
and  dynamic  variables  we  can  prescribe  the  boundary  conditions  of  the 
problem  as  follows 

9F,  9F, 

v^-^  +  w-^-  =  0  on  Fl(y,  z)  =  0  ,  (10) 

and 


p  =  p  =  const.  on  F?(y,  z)  =  0  ,  (11) 

where  F^(y,  z)  =  0  denotes  the  wetted  body  surface  and  F^(y»  z)  =  0, 
the  free  streamlines  that  form  the  wake  boundary,  but  are  unknown 
a  priori.  The  constant  pressure  Pc  inside  the  wake,  called  the  wake 
underpressure  (pc  <  p^  ,  the  pressure  at  infinity,  in  general)*  is 
assumed  to  be  given.  In  order  to  make  possible  our  consideration  of 
the  general  case  of  pc  <  p^  ,  a  theoretical  model  of  finite  wake  is 
required,  for  which  we  again  adopt  the  'open-wake  model'  as  in  the 
previous  case  (of  uniform  free  stream  treated  in  Report  AMT-R-101). 
Finally,  we  have  the  following  boundary  conditions  at  infinity: 


U, 


w 


0, 


00 


as 


y2  + 


QC 


(12] 


3.  1  Small  perturbation  theory 

To  render  this  nonlinear  problem  tractable  to  complex  functional 
analysis,  we  shall  assume  that  the  primary  flow  is  only  weakly  sheared, 
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V!*  -~v  ^ .-*-  x 


i.  e.  ,  •y  <  «  1  as  measured  with  respect  to  the  body  dimension  i  .  Thus 

we  seek  the  solution  in  the  form  of  an  asymptotic  expansion  for  small  yi, 
or  in  complex  notation 
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x  =  y  +  iz  ,  x  =  y  -  iz  ,  (13 ) 

W(x,x;y)=  V  -  iw  =  Wo(x)  +  Uyz  +  Y  Wj(x)  +  0(y2)  , 

(14) 

Wn(x)  =  dfn(x)/dx  =  vn  '  i  wn  .  (n=0,  1,2,  ..)  (15) 

where  fn(x)  =  $n(y>  z)  +  l|4,n(y*  z)  (n  =  0,  1,  2,  .  .  .  )  are  analytic 

functions  of  x»  regular  everywhere  in  the  interior  of  the  flow  field,  but 
may  have  branch-point  type  singularities  at  isolated  points  on  the  flow 
boundary.  The  wetted  body  surface  will  be  given  parametrically  by 

y  =  y(s),  z  =  z(s)  (0  <  s  <  S)  ,  (16) 

where  s  is  the  arc  length  along  the  body  surface,  running  from  point  A 
(where  s  =  0)  to  point  B  (where  s  =  S,  see  Fig.  1).  To  facilitate  our 
functional  analysis,  we  continue  to  employ  the  same  parametric 
4  =  £  +  in  plane  as  in  the  case  of  uniform  current,  such  that  the  flow 
region  is  mapped  conformally  into  the  lower  half  of  the  L, -plane,  with  the 
wetted  body  surface  mapped  onto  (-1  <  £  <  1,  r)  =  0),  the  stagnation 

point  provisionally  onto  the  origin  t,  =  0,  and  the  cavity  boundary  onto 
(|£|  >  1,  r|  =  0).  Substituting  the  expansion  (14)  in  (10)  and  using 
(3Fj /dy)/ (dF^/dz)  =  -z' (s )/y' (s ),  we  obtain 

(vo  +  YUz  +  VVj  +  ...  )z'  (s)  =  (wo  +  v  Wj  +  .  .  .  )y'(s)  ( |£  i  <  1,  ii=0), 


tv 
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The  first  order  theory  has  been  developed  to  provide  solutions  as 
already  presented  in  AM  Tech  Rept  AMT-R-101.  For  the  convenience  of 
applications  in  the  sequel  we  delineate  the  main  steps  of  the  analysis  as 
follows.  First,  we  represent  the  finite  wake  by  the  generalized  open-wide 
model  introduced  by  Wu  &  Wang  (1964a),  according  to  which  the  conformal 
transformation  between  the  leading  order  complex  potential  fQ  =  +  i4>o 

and  the  parametric  variable  £  =  £  +  iq  (as  depicted  in  Fig.  2)  is  given  by 
the  Schwarz-Christoffel  mapping. 


df 

o 


U4-£c> 


{S-S0)<S-S0> 


(21) 


where  A  is  a  real  constant,  £o  is  the  image  point  of  ^  =  oc  (also  fQ  =  oo), 

its  complex  conjugate,  £  is  a  real  constant  marking  the  rear  end  of 
o  c 

the  near  wake  at  which  the  constant  wake  pressure  region  terminates. 

To  best  effect  satisfaction  of  the  kinematic  and  dynamic  conditions 
we  introduce  the  logarithmic  hodograph  variable 


T 

O 


qc  being  the  constant  flow  speed,  to  leading  order,  on  the  wake  boundary 
Here  6q  is  the  leading-order  angle  subtended  by  the  flow  velocity  to  the 
y-axis.  In  terms  of  frQ(4)  =  T0(4»Tl)  +  i0o(|,ri)  conditions  (17a)  and  (19a) 
become 


0~o(t)  E0o(£.O')=  *#(-!)  +  P(s<£)) 


(IS  I  <  1) 


(23) 


0 


del  >  i).  (24) 


V*>  £  To^'°  )  = 

wherej=f  (-£,)  denotes  Heaviside's  step  function 

#(-£)  =1  for  |  <  0  ,  #(-|)  =  0  for  4>0, 

and 

p(s)  =  tan-1  [z*(s)/y'(s)]  .  (25) 

This  Riemann- Hilbert  problem  can  be  resolved  as  follows.  We  first 
continue  analytically  the  flow  field  in  the  lower  half  4-plane  into  the  upper 
half  4 -plane  by  the  Schwarz  reflection 

syr)  = -"Tyrr .  m 

where  the  overhead  bar  denotes  taking  the  complex  conjugate,  or 
equivalently, 

to(4,-»|)  =  -to(£,ti)  ,  eo(4,-q)  =  eo(4,T!)  .  (27) 

Then,  noting  that 

«0*(£>  =  7t^,  +  ie±(S>  =  +  To  (£>  + ieo(S)  ’  (28) 

in  which  the  last  step  follows  from  (2  7),  conditions  (23)  and  (24)  can  be 
written  as 

fio  +  fio  =  (given  for  |£|  <  1)  ,  (29) 

U+o  '  %  =  2To{e)  =0  (lei  >  1)  -  (30) 

The  associated  homogeneous  problem  (corresponding  to  (29)  and  (30), 
is  defined  by 


a 

\ 

% 

* 

% 

J 

> 

f 

\ 
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s. 

Ht(£)  +  H'(£)  =  0 


Ht(S)  -  H'(£)  =  o 


del  <  i)  . 
del  >  i)  > 


where  H(4)  is  required  to  be  holomorphic  everywhere  in  the  £, -plane 
except  possibly  at  4  =  ±  1  and  4  =°°i  it  has  infinitely  many  solutions: 


2  1/2 

H(4)  =  (4  -  1)  '  R(£)  » 


where  R(4)  is  any  suitable  rational  function  of  4-  Here  the  function 


2  1/2 

(4-1)  is  defined  with  a  branch  cut  in  the  4-plane  from  4  =  -1  to  1 


along  r|  =  0  so  that  it  is  single- valued  in  the  entire  cut  plane  of  4  and 
tends  to  (,  as  |  4  |  oo  for  all  arg  4  •  We  next  introduce  the  new  function 


F(4)  =  fco(4)/H(4)  - 


From  (2 9 ) - (3 2 )  we  see  that  F(4)  has  a  jump  across  the  £-axis  by 


F+(£)-F"(£)  =  (l£  +  flj/H  +(£)  =  2i6^  (£)/H+(£)  (|£|<1), 

(35) 


=  (<  -  «j/H+(£)  =  0  (lei  >  1)  . 


o  o 


A  particular  solution  can  therefore  be  obtained  for  F(4)  by  applying 
Plemelj's  formula  (see,  e.  g.  Muskhelishvile  1953,  §17) 


=  _J_  f°°  Flm.: 

2-rri  J  t  - 


from  which  we  deduce 


a!-!'-* 


I.  •  .  1  .  -  •  ■- 3_ ■».  *1  ..•*  V *  «. 


»(£)  =F{4)H(£)  =  -H(4)\  — -  dt  .  (38a) 

°  -i  «(*>(*-&> 

Under  the  required  condition  that  ft0(C)  he  analytic  and  regular  everywhere 
in  the  cut  4-plane,  in  particular  at  4  =  +_  1  and  4  =  00  *  the  only  possible 
form  for  H(4)  is 

H(4)  =  (42  -  1)1/2  ,  (38b) 


and  then  (38a)  gives  the  unique  solution  of  fiQ(4)*  Using  the  boundary 
value  (23)  for  8 ~  (t)  in  (38)  and  carrying  out  some  integration  finally  yields 


n(C) 


log 


i+i(42-i)1/2 


TTl  J 

-1 


21/2 

£-4> 

i-t 


P(t)dt 


(39) 


where 


P(S)  =  p(S(e» 


del  <  i)  .  (40) 


The  solution  is  thus  determined  to  leading  order  in  parametric 
representation,  £q  =  f q(4 )  ,  ftQ  =  fto(4),  given  by  (21 )  and  (39).  The 
physical  plane  is  then  given  to  leading  order  by  quadrature 

zo(U  =  I  vrtT)  ZT"d^  (41) 

The  arc  length  sq(4),  measured  along  the  body  surface  from  |  =  -1,  is 


f 


rc(e) 


v(|)d^ 


Ml) 


A 


(-i  <  e  <  i)  , 


(42) 


W£)  =  [1+  (1-0  ]  - 2 - zr~2 -  (44) 

(€-40)  (l-&0) 

The  integral  in  Eq.  (43)  assumes  its  Cauchy  principal  value.  The  total 
wetted  arc  length  is  specified  by 

S  =  so(l)  .  (45) 


Moreover,  the  velocity  condition  at  x  =  oo  ,  or  t,  =  requires  that 


a0K0>  =  log  ir=  ilog  (1  +  CT)  • 


(46) 


This  completes  the  formal  solution  to  leading  order  for  the  wake  flow 
past  a  ship  section  with  a  uniform  incident  current.  The  foregoing  analysis 
provides  a  set  of  nonlinear  functional  equations  (42)  and  (46)  which  we  write 
symbolically  as 


so(£)/(A/qc)  =  £  [so(|),  (3(so(£));  t,Q] 


where 


to  =  (u/qc)?[*o(i). 

&  [s0(5).  P  (•„(«));  tj 


P  C0(6))i  tj  . 


C 

■r 


r  (£> 


v(4;40)de  . 


(47) 

(48) 


e 


(49) 


o 


3[80(€).P(a0(g));&0]  =  [i+iao-D  ]exP^j  <— 2>  — — 

1  1  t-L 


Solution  of  the  above  set  of  nonlinear  functional  equations  can  be  obtained 
by  numerical  methods.  Further  development  of  a  workable  numerical 
method  notwithstanding,  we  proceed  to  discuss  below  the  first-order 
theory  by  which  the  effect  due  to  a  weak  vertical  shear  of  the  incident 
current  can  be  evaluated. 

3.  2  The  first-order  perturbation  theory  for  a  weakly  sheared 


current 

To  first  order  in  the  shear  gradient  y,  condition  (17b)  -  after 
eliminating  z  '/y'  in  it  by  making  use  of  (17a)  -  and  condition  (19b) 
become 


v  w,  -  w  v,  =  Uw  z 
o  1  o  1  o 


(l£ !  <  q  =  °") , 


v  v,  +  w  w,  =  U(Uz  -  v  z)  (  |  >1,  q  =  0  ) 

ol  ol  '  oo  o  '  1  1  ' 


These  two  equations  can  be  recast  as 


Im  G, 


Uw  z 
o  o 


(  i£  I  <  -9  =  0)  , 


Re  G,  =  U(Uz  -  v  z  ) 
1  '  oo  o  o' 


(lei  >  1,  q  =  o’)  , 


where 


Gl(e)  =  =  (vo  +  iwo)(Vj  -  iWj) 


In  (51)  and  (52)  we  have  replaced  the  unknown  function  z(e,q)  by  its 
leading-order  solution  (41).  Now  with  the  function  Gj(£,)  analytically 
continued  into  the  upper  half  t, -plane  by 
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(62) 


We  further  note  that  the  stagnation  point  on  the  body  surface  is  no  longer 
necessarily  at  £  =  0  but  is  generally  shifted  to  l,  -  on  the  real 

£,-axis  due  to  the  contribution  from  Wj,  i.  e.  , 

W0(£*>  +  =  0  ,  (63) 

being  clearly  of  order  y  .  This  in  turn  requires  that  the  transforma¬ 
tion  (21)  be  modified  to  the  following  form 

df(i)  (4  -  4*><4  -  !c> 

- -  =  A  - 2 - 3-2  ,  (64) 

<4  -  40)  (4  -  40) 

in  which  the  relocated  zero  of  df^Vd4  at  t,  =  £j!c  is  required  to  make  the 
new  stagnation  point  a  branch  point  of  the  stagnation  streamline.  With 
this  up-graded  f^(4)»  the  physical  plane  can  be  calculated  anew  by 
integrating 

df(1)/dx  =  Wo(4)  +  vWjd)  ,  (65) 

giving 

(1 )  r  &  1  df^  ^ 

X  =  J  W  (4)  +  YWj(4)  "dT-  (66) 

Finally,  the  same  conditions  as  (42)  and  (46)  for  the  leading  order  solution 
can  be  applied  to  ^  of  (66)  and  ^  of  (62);  and  the  remaining  task  is  to 


develop  a  numerical  method  for  computing  the  solution. 


16 


u 


4.  The  hydrodynamic  forces  and  moments 

The  hydrodynamic  forces  and  moments  of  significance  during  the 

motion  of  a  moored  ship  are  primarily  the  sway  force  F^  and  yaw  moment 

M  .  To  facilitate  possibly  accounting  for  large  excursions  in  the  horizontal 
z 

plane  we  employ  a  coordinate  system  (x,  y,  z)  fixed  with  respect  to  the  ship. 
The  remaining  three  modes  of  rigid-body  motion,  namely  the  heave,  pitch 
and  roll  can  be  assumed  to  be  unimportant  in  moored  vessel  problems,  as 
is  generally  the  case,  and  are  accordingly  ignored. 

The  motion  of  a  moored  ship,  even  for  the  case  of  a  single-point 
mooring,  is  relatively  slow,  as  compared  to  vertical  ship  motions  in  waves, 
due  to  the  predominant  inertial  effects.  Thus  it  is  reasonable  to  neglect  the 
wave  effects  associated  with  the  wave  making  by  the  current,  which  can  be 
evaluated  separately  later.  For  the  same  reason  we  assume  that  the 
transverse  flow  with  sway  and  yaw  can  be  analyzed  on  the  basis  of  the  low- 
frequency  limit  of  the  free- surface  conditions.  In  linear  theory  the  two 
free-surface  conditions  (when  expressed  with  respect  to  a  coordinate 
system  moving  with  the  current)  for  small  free  surface  elevation  q(x,  y,  t) 
and  velocity  potential  ^(x,  y,  z,  t)  read 


+ 

at 

<JQ 

-3 

II 

o 

(2=0)  , 

(67) 

M  - 

3z 

9ri 

at 

(2=0)  . 

(68) 

For  time  harmonic 

motions. 

with  the  factor  elw*,  the 

above  two  conditions 

combine  to  give 

d4> 

dz 

2 

= 

g 

0  .  (z=0) 

(69) 
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In  the  low-frequency  limit,  oo  i/g 


0,  the  above  condition  reduces  to 


9£ 

3z 


0 


(z  =  0)  ,  (70) 


that  is,  the  free  surface  becomes  in  this  limit  a  rigid  horizontal  plane. 

With  this  approximation  the  flow  field  can  be  continued  into  the  upper  half 
space  z  >  0  by  mirror  image,  and  the  problem  becomes  one  of  irrotational 
flow  of  an  inviscid  fluid  past  a  rigid  double  body,  symmetrical  with  respect 
to  the  plane  z  =  0.  Thus  the  wetted  portion  of  the  hull  is  merely  the  low*»r 
half  of  the  double  body. 

In  our  analysis  of  the  transverse  hydrodynamic  force  associated 
with  the  total  transverse  current  velocity 


V  (x,  z,  t)  =  Uc(x,t)(l+yoz)  -  UG(t)  -x  f2(t)  =  Uyc-xS,  (71) 

which  now  includes  slowly  varying  spatial  and  temporal  variations,  we 
distinguish  between  two  fundamentally  different  flow  regimes.  One  is 
characterized  by  a  viscous  wake  formation  that  occurs  at  sufficiently  high 
incidence  angles  of  the  current;  the  other  takes  place  at  low  incidence  angles 
when  the  lateral  wake  swings  backwards  to  merge  with  any  existing  vortex 
sheet  shed  from  the  keel  into  a  trailing  stern  wake.  For  practical  purposes 
a  pragmatic  continuous  transition  between  the  two  regimes  can  be  sought  by 
numerical  matching. 

For  the  case  of  lateral  wake  flow  problem  we  first  evaluate  for  the 
local  strip  (at  station  x  and  at  time  t)  the  two-dimensional  wake  flow  in  a 
vertically  sheared  current  following  the  analysis  and  solution  method 
presented  in  the  previous  section.  The  total  y-component  velocity  given 
by  (71)  can  be  recombined  as 
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V(x,z,t)  =  U2(*t)  [1  +  Y  Z1  >  (72) 

in  which  the  specific  shear  rate  y  maY  be  taken  as  general  as  a  slowly 
varying  function  of  x  and  t.  By  integrating  the  resulting  pressure  over 
the  wetted  hull  of  the  double  body  with  respect  to  the  wake  underpressure 
as  datum,  we  obtain  the  y-component  of  the  hydrodynamic  differential 
force  as 

dFyc  =  }pwCDC(x.t;or)U2(x,t)h(x)dx  ,  (73) 

where  h(x)  is  the  sectional  draft  of  the  ship  (see  Fig.  3).  This  result 

is  now  applicable  to  cases  more  general  than  (17)  of  Part  I,  since  both 

spatial  and  temporal  variations  in  the  current  velocity  have  been  taken 

into  account,  provided  that  these  variations  are  slowly  varying  and  the 

water  depth  is  not  much  shallower  than  twice  the  ship's  draft.  The  total 

transverse  current  force  and  the  yaw  moment  can  then  be  obtained  by 

integration  from  the  ship's  bow  end  at  x  =  -I,  to  the  stern  at  x  =  f  as 

1  b  s 

s 

Fyc(t)  =IPwJ  CDC(x,t;a)U^(x.t)h(x)dx  ,  (74) 

-*b 

Mc(t)  =  I  pw  y  CDC  (x,t;a)U*(x,t)h(x)xdx  .  (75) 

Wo  remark  again  that  the  end  strips  at  the  stern  end,  x  =  f  ,  and  at 
the  bow  end,  x  -  can  ^an<^^e<^  by  employing  the  refined  scheme 

described  below  Eq.  (17),  Part  I. 

For  the  case  of  small  incidences  of  the  incoming  current  with 
respect  to  the  ship's  centerplane,  the  general  theory  of  Wu  (1971) 


HA 


developed  for  unsteady  slender-body  flow  remains  valid  even  when  the 
current  has  spatial  and  temporal  variations  in  velocity  provided  they  are 
slowly  varying  in  units  of  the  body's  transverse  dimension  and  its 
characteristic  time  of  motion.  Thus  the  local  transverse  force  (called 
the  differential  lift  due  to  circulation)  acting  on  a  body  strip  of  unit 
thickness  (see  Wu  1971)  is  given  by  <£nir.  =  dL  .  /dx,  where 


=  +  Ui  t)m22^x^  (without  side  vortex)  ,  (76a) 


m22(x)(9^  +  ui  hm*.*) 


'1  9x' 


(with  side  vortex)  ,  (76b) 


where  Uj(x,  t)  is  the  local  longitudinal  velocity  component  with  respect 
to  the  body  section  and  m22(x)  is  the  sectional,  frequency- independent 
added-mass  coefficient  of  the  lower  half  of  the  double  body  for  the  sway  mode. 
Strictly  speaking,  the  above  formulas  for  the  differential  lift  holds  only 
for  currents  having  no  vertical  gradient  of  shear,  that  is,  with  y  =  0 
in  (72).  AMien  y  ^  0,  attempts  may  be  made  to  derive  an  equivalent 
added-mass  coefficient  to  replace  m^ix)  in  (76),  but  such  a 

correction  for  the  vertical  shear  effect  would  require  considerable  amount 
of  new  developments,  analysis  and  calculation. 

On  the  same  premise  the  differential  induced  drag  associated  with 
the  local  lift  of  this  circulatory  flow  can  be  obtained  by  applying  the  general 
slender-body  theory  of  Wu  (1971)  to  yield 

Air  =  dDcir/d*  "  U1  -  *2n2>m22(x)]  .  <77a) 


1  /  W  9  ^  9  WTT2  2cA 

Im22 <X>(9F  +  U1  8^)(Uvc  -  x  n  >  ’ 


(77b) 


according  as  there  is,  (77a),  or  there  is  not,  (77b),  a  local  vortex  sheet 
being  shed  from  the  keel  of  the  hull  surface.  The  total  longitudinal  and 
transverse  (i.  e.  ,  the  x-  and  y-component)  hydrodynamic  forces  are 
therefore  given  by  the  integrals 


(78) 


(79) 


where  the  local  incidence  angle  Cl  may  also  be  a  slowly  varying  function 
of  x  and  t.  Finally,  we  can  express  the  yaw  moment  about  the  origin 


:„(x.  t)  x  dx 


5.  The  rudder  force 

The  general  slender-body  theory  can  further  be  applied  to  predict 
the  rudder  force  on  the  assumption  that  the  hull-rudder  system  as  re¬ 
presented  by  its  double  model  is  a  single  low-aspect-ratio  lifting  surface 
with  prescribed  surface  displacement.  As  the  contribution  from  the 
transverse  current  velocities  to  sway  force  and  yaw  moment  have  already 
been  determined  for  both  cases  of  finite  and  small  values  of  the  incidence 

angle,  we  need  onlyconsider  the  action  of  the  rudder  as  a  deflection  of  the 
stern  part  of  the  slender  hull,  which  is  the  rudder  itself  and  its  double 
model  image,  from  its  stretched-straight  position  by  the  amount 


y  =  -<=k<t><*  -  JRI  «R  <  x  <  /s.  h,R<  U|  <  h2R) 

=  0  (elsewhere)  (81) 

where  the  above  specified  region  in  x  and  |z|  denotes  the  rudder  plan- 
form,  x  =  i R  is  the  coordinate  of  the  rudder  axis  (assumed  vertical)  and 
(a^(t)  is  the  angle  of  rudder  deflection.  This  rudder  displacement  gives 
rise  to  an  equivalent  transverse  flow  velocity  (relative  to  body) 

V(x,  t)  =  (gf  +  Uj  g^)@R(t)(x  -  fR)  =  +  feR(t)(x  -  fR)  + 

(82) 

Substituting  this  relation  for  V  in  (76)  and  integrating  the  resulting  lift 
over  the  rudder  section  gives  the  rudder  force 

FyR  =  Z  ^RCR  +  2UlCR^k  +  U1  ^R^  mR  '  (83  > 

where  CD  =  i  -  i_  is  the  constant  chord  of  the  rudder  and  m_  is  the 

xv  S  xv  K. 

two-dimensional  added  mass  of  the  rudder  which  is  assumed  to  be  a  con¬ 
stant  along  its  chord.  The  yaw  moment  due  to  rudder  action  is  simply 

MR  =  FyRiR  ’  (84) 


This  contribution  from  rudder  to  the  sway  force  and  yaw  moment  can  be 
combined  with  that  from  other  sources  in  assembling  the  equations  of 
motion  for  the  moored  ship. 
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6.  Wave  radiation  from  a  ship  moored  in  steady  current 

The  problem  of  interaction  between  the  ship  hull  and  a  steady 

current  has  been  resolved  in  the  preceding  sections  for  the  low  frequency 

limit,  that  is,  with  the  assumption  that  the  gravity  effect  on  wave  making 

by  the  ship  held  fixed  in  the  current  is  of  secondary  importance  and  was 

therefore  neglected  in  order  to  facilitate  a  full  account  of  the  primary 

effect  of  wake  formation.  (To  give  a  rough  estimate  for  comparison,  we 

note  that  for  the  Froude  number  of  a  typical  ship  of  length  l  to  be 

1  /2 

Fr  =  U/(gi)  <  0.2,  the  wave  resistance  coefficient  is  <  0.  1  as  a 
rule  whereas  the  pressure  drag  coefficient  associated  with  a  transverse 
wake  due  to  a  cross  current  can  be  as  large  as  1  or  somewhat  greater.  ) 

We  now  proceed  to  investigate  this  gravity  effect  and  give  a 
separate  account  of  the  energy  loss  due  to  wave  generation  caused  by 
interaction  between  the  hull  and  a  steady  current.  We  shall  assume  that 
the  incoming  stream  is  a  steady,  at  least  a  quasi- steady  current,  without 
carrying  any  incident  gravity  wave,  and  that  the  incidence  angle  Ot 
(between  the  current  velocity  and  the  longitudinal  x-axis  of  the  ship)  is 
arbitrary  (0  <  ot  <  360°).  The  last  assumption  distinguishes  the 
present  problem  from  the  classical  theory  of  wave-making  resistance  of 
ships;  in  fact,  we  note  that  the  vast  literature  is  very  much  limited  to  the 
case  of  Ot  =  0  as  only  this  case  is  required  for  consideration  of  ship 
propulsion. 

To  effect  a  viable  engineering  solution  to  the  present  problem, 
we  shall  assume  that  when  the  transverse  wake  is  manifest,  the  shape  of 
which  is  supposed  to  be  already  determined  by  the  method  proposed,  the 
wake  boundary  plus  the  front  part  of  the  wetted  hull  outside  the  wake  will 


be  taken  as  the  new  'floating -body'  surface  that  is  responsible  for  wave 
generation.  This  assumption  is  made  on  physical  ground  that  the  wake 
boundary  (assumed  stationary  in  the  mean)  is  a  stream  surface,  but  its 
enforcement  will  make  the  new  body  no  longer  as  thin  or  slender  as  the 
original  ship  hull.  Consequently  the  classical  thin-ship  theory  and 
slender-ship  theory  should  both  become  nonapplicable  for  the  present 
case,  but  the  method  of  Green's  function  can  still  be  useful.  To  render 
the  problem  tractable,  however,  we  shall  further  assume  that  the  resulting 
waves  are  small  in  amplitude  so  that  the  free- surface  boundary  conditions 
can  be  linearized.  Thus  we  have  the  following  formulation. 

There  exists  an  additional  velocity  potential  #(x,  y,  z)  which  is 
solely  responsible  for  wave  generation.  This  $(x)  must  satisfy,  in  the 
linearized  flow  region  0  >  z  >  -hQ  =  const.,  the  following  equations 


6  +0  +  <j>  =0 

xx  yy  zz 

(-h  <  z  <  0)  , 

l 

(85) 

<t>  +  K  <t>  =  0 

XX  z 

(z  =  0,  K  =  g/U*  )  , 

(86) 

=  -U  n. 

3n  cl 

(x  on  SBW)  , 

(87) 

2  2 

<t>  =  o((x  +  y  )  ) 

2  2 

asx+y  -•■oo,  x<0. 

(88) 

6  -  0 
^z 

(z  -  -v  • 

(89) 

Here,  we  have  adopted  a  new  rectangular  coordinate  system  with  the  x-axis 
pointing  in  the  direction  of  the  current  velocity  U^,  the  z-axis  remaining 
the  same  as  before  so  that  the  z  =  0  plane  coincides  with  the  undisturbed 
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water  surface  S^,  the  z  =  -1^  plane  is  the  horizontal  floor  of  the  water 
layer  (including  hQ  =  oo  as  the  deep  water  limit),  g  is  the  gravitational 
acceleration.  Equation  (86)  is  the  linearized  free-surface  condition 
which  is  the  stationary  version  of  (69).  In  (87),  d<t>/ 9n  =  n  •  ,  and 

n^  is  the  x-component  of  the  outward  (from  the  fluid)  normal  n  to  the 
surface  S„,.r  =  +  S,.r  which  is  the  union  of  the  wetted  body  and  wake 

boundary.  The  order  limitation  (88)  is  the  Sommerfeld  radiation  condition 
for  three-dimensional  water  waves.  The  free  surface  elevation  Z(x,  y)  is 
determined  by 

U 

Z(x,  y)  =  — C(^x(x,  y,  0)  .  (90) 

The  hydrodynamic  force  acting  on  the  hull  due  to  the  effects  of  wave 
making  is 

F  =  -puj’  0xndS  +  (pw-pc)ndS  ,  (91) 

SB  SW 

where  n  is  the  unit  normal  pointing  into  the  hull  and  the  wake,  and 
p  =  pw  on  the  wake  boundary  enclosing  the  near  wake  with  a  specified 
base  pressure  pc> 

The  desired  Green's  function  may  assume  the  form 

G(x;£)  -  r"1  +  G*(x;£)  ,  (92) 

whc  r  e 

2  2  2  1/2 

r  =  -  £  j  =  [(x  -  £)  +  (y  -  tj)  +  (z  -  £,)  ]  (93) 

is  the  distance  between  two  field  points  x  and  £  ,  G  is  harmonic  and 
regular  in  the  flow  region  0  >  z  >  -ho, 


V  G 


0 


(94) 


■r  -r  --  w  /.v  .  .  v  r 
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(0  >  z  >  -h  )  . 

'  o 

The  function  G  satisfies  the  following  boundary  conditions 
G^(x;e.ri.O)  +  KG^fx  ;4,  tj.  0)  =  0 
G^feye.^-h)  =  0 

and 

-1/2 

G  =  0(  +  x\“  |  )  as  4^  +  r\“  -*  oo  ,  |<0. 


(95) 

(96) 

(97) 


We  now  apply  Green's  third  formula  to  the  region  bounded  by  the 

plane  S ^(t,  =  0),  the  bottom  S^(4  =  -1^)»  a  large  control  surface  2, 

2  2  1/2 

(4  +  ri  )  =  R,  the  wetted  hull  surface  Sg  and  the  wake  boundary  S^. 

Then  we  have 

4Tr^(x)=^[|f  G(x;|)  -  <p(t)  |^]  dSg  (98) 

S 

where  S  =  S£  +  +2+  Sg  +  Sw  ,  and  dS|  signifies  a  surface  element 

at  |  on  S.  By  assuming  for  $  and  G  the  following  asymptotic 
representations  (taking  the  real  component  for  physical  interpretation) 

_1  e+iKrf(0)  +  0(i}  as  r  _  00<  4  >0  (99) 

)  r 

where  1(6)  stands  for  a  function  of  6  =  arctan  (•»-)/ 4 )  only,  as  can  be 


verified  a  posteriori,  it  is  readily  shown  that  the  integral  over  the 
control  surface  2  vanishes  as  R  -►  oo  .  After  using  the  boundary 
conditions  for  <j>  and  G,  we  obtain 


4ir*<x)  =  j  C-Uixj  (£)G(x;| )-  *(£)  J  <*S  +  ^  j  ~  (</>G^  -  <^G)d£dT, 


SD  +  S 
B  w 


The  last  integral  over  the  free -surface  may  be  intregrated  by  parts 


4ir0(x)  =  [Unj  (£  )G(*;£ )+  <p  (£ )  ff  ]dS  +  i  <£  G  -  <p  G^  )dr,  , 


S„+S 
B  w 


(100) 


where  the  line  integral  is  taken  counterclockwise  around  the  contour 
in  the  £-tj  plane  which  is  the  intersection  of  Sg  and  S£.  (The  line 
integral  around  the  intersection  of  and  2  vanishes  as  R  -*•  oo  .  )  Now 
if  in  the  surface  integral  we  let  a  field  point  x  -*•  x  Q  on  Sg+  Sw„  then 
by  a  well  known  theorem  (a  lemma  of  Green's  theorems)  the  integral 
with  8G/9n  converges  to 

xU™x  I  *(l>Gn&ii)dSe  =  '2^(^0)+^  ^(l)Gn^o;l)dS5  • 


We  then  have  the  following  integral  equation  for  <p(x),  where  x  is  a  point 
of  SD  : 


2ir*>(x)+C  ip(  |)G„(x;£)dSc  =  -U  \  n,  (|  )G(x  ;£  )dSt  +  4  <£  (0,  G-^GC  )dq. 


by  numerical 


When  this  equation  can  be  solved  for  <f>  for  x  on  Sgw 
methods,  then  <j>  is  determined  everywhere  by  (100).  However,  for  the 
general  case  of  a  blunt  body  in  transverse  current,  no  numerical  method 
is  known  to  exist  and  it  may  require  development.  The  Green's  function 
satisfying  (94)-(97)  for  infinite  depth  (hQ  =  oo)  can  be  given  as  follows: 


G(x;£)  = 


-1  -1 
r  -  rx 


(102) 


tt/2 


sec^0d0 


oo 

§ 


,k(z  +  £,) 


k-Ksec  0 


2 —  cos  k^(x-|)cos  ^(y-qjdk 


tt/2  2 

-4K^  eK(z+4)sec  ^sec^Q  sin[Kj  (x-|  )sec20]cos[Kj  (y-q )seC20]d0 


where 

1  /2 

=  [(x-^)2  +  (y-t])2  +  (zH)2] 


(102b) 


kj  =  k  cos  0,  k2  =  ksin0,  Kj  =  Kcos  0,  K2  =  Ksin  0.  (102c) 


7.  The  shallow  water  effects 

In  investigating  the  shallow  water  effects  on  the  hydrodynamic 
characteristics  of  a  ship  moored  in  an  incident  current  we  consider  the 
possible  resulting  changes  in 

(i)  the  transverse  wake  formation  and  the  associated  pressure 
drag  (which  is  a  quasi-steady  force); 

(ii)  the  added  mass  coefficients  for  (unsteady)  sway  and  yaw 

modes  of  motion,  especially  in  the  presence  of  a  transverse 
wake; 


(iii)  the  energy  loss  due  to  wave  radiation  and  the  resistance  due 
to  wave  making. 

Before  we  proceed  with  the  formulation  of  these  problems,  it  is  worth 
noting  that  the  main  problem  at  hand  distinguishes  itself  from  the  more 
classical  problem  of  propulsion  and  maneuvering  of  a  ship  in  shallow 
water  by  having  a  transverse  wake  as  a  drastic  new  feature  for  a  ship 
moored  in  current.  The  nonlinear  effects  that  arise  from  the  transverse 
wake  formation  now  combine  with  the  nonlinear  shallow  water  effects  to 
make  the  problem  extremely  difficult  and  there  is  practically  no  literature 
available  that  are  related  to  the  problem  we  now  seek  to  resolve.  The 
complexity  is  somewhat  alleviated  by  the  operating  condition  that  the 
depth  Froude  number  of  the  primary  current  usually  lies  in  a  low  sub- 
critical  range  (less  than  0.  5  say).  Consequently,  the  nonlinear  boundary 
conditions  on  the  water  surface  can  be  linearized,  and  this  will  be  the 
basis  of  our  subsequent  analysis  of  the  problem. 

7.  1  The  transverse  wake  and  its  pressure  drag 

Strip  theory  can  be  applied  to  study  of  the  shallow  water  effects  on 
the  transverse  wake  formation  and  its  corresponding  pressure  drag.  We 
shall  however  assume,  as  explained  above,  that  the  Froude  number  based 
on  the  water  depth,  F^,  of  the  primary  current  of  velocity  Uc,  lies  in  a 
low  subcritical  range,  i.  e.  , 

Fh  =  Uc/(gho)1/2  <  0>5  '  (103) 

where  hQ  is  the  (uniform)  water  depth.  Adopting  strip  theory  for 
evaluating  the  local  two-dimensional  flow,  as  in  section  3,  and  restricting 
our  attention  to  the  low-frequency  limit,  the  free -surface  condition  can  be 
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approximated  by  the  rigid-surface  condition  80/9z  =0  on  z  =  0,  while 
the  same  condition  also  applies  on  the  bottom  x  =  -hQ.  Thus  the  wake 
flow  becomes  equivalent  to  one  past  the  double  model  placed  symmetrically 
in  a  closed  tunnel  of  spacing  2hQ  — which  is  a  pure-drag  problem. 

The  wall  effects  in  the  pure-drag  case  of  plane  wake  flow  past  an 
arbitrary  body  held  in  a  closed  tunnel  have  been  investigated  by  Wu  et  al. 
(1971),  both  theoretically  and  experimentally.  It  has  been  shown  there 
that  for  wake  flows  past  fairly  blunt  bodies,  such  as  the  double  model  section 
of  a  typical  mid- ship  hull,  the  wall  effects  are  generally  very  small  and 
can  be  neglected  in  predicting  the  pressure  drag  coefficient  if  the 

comparison  is  based  on  the  same  wake  underpressure  coefficient  a.  For 
example,  when  the  water  depth  is  no  less  than  1.  5  times  the  ship's  draft, 
the  drag  coefficient  will  have  only  a  small  reduction  from  its  deep¬ 

water  limit,  by  a  fraction  no  greater  than  3  %for  a  up  to  as  large  as  2. 

Such  predictions  have  been  strongly  confirmed  by  experiments.  Since  a 
will  be  employed  in  the  present  model  as  an  empirical  parameter,  we 
may  therefore  assume  the  shallow  water  effects  on  the  transverse  wake 
pressure  drag  to  be  negligible. 

7.2  The  added-mass  coefficients 


As  already  explained,  the  primary  effects  of  shallow  water  on  the 
unsteady  hydrodynamic  characteristics  are  reflected  in  the  resulting 
added-mass  coefficients  for  sway  and  yaw  motions.  The  leading-order 
results  for  these  motions  can  be  constructed  from  a  low-frequency 
strip-theory  analysis  such  that  the  sectional  two-dimensional  added-mass 
coefficient  is  evaluated  for  the  double  model  placed  symmetrically 

in  a  channel  of  total  depth  2h  .  When  the  transverse  flow  has  no  wake 


(that  is,  when  flow  separation  is  negligible)  m^fc)  can  be  determined 
by  classical  methods,  such  as  the  Green's  function  method,  the  boundary 
integral  equation  method  or  conformal  mappings.  For  given  hull  section 
profiles,  can  be  evaluated  by  these  numerical  methods.  The 

general  trend  of  increase  in  rn£2  ^ue  bottom  confinement  can  be 

exemplified  by  the  flow  past  a  circle  of  radius  a  placed  midway  between 
two  plane  walls  at  distance  2hQ  apart,  for  which  we  can  show  by 
perturbation  theory  that 

?  2  2 

m22  “  p7ra  t1  +  12  ~~2  +  J  »  (104a) 

ho 

o 

for  a/hQ  «  1.  If  the  circle  is  replaced  by  a  flat  plate  of  breadth  2a,  held 
broadwise  to  the  bounded  stream,  we  have 

?  2  2  . 

m22  =  p7ra  ^  ^  U  ^  2  +  ]  •  (104b) 

^ ‘  24  h  o 

o 

These  results  show  that  the  sway  and  yaw  added-mass  coefficients  are 
both  increased  by  proximity  to  the  bottom  (i.  e.  ,  as  water  depth  is 
decreased)  and  that  the  increment  seems  to  depend  more  sensitively  on 
hull  shape  than  the  leading -order  quantity. 

The  above  approach,  however,  is  further  complicated  in  practice 
by  the  presence  of  a  transverse  wake.  On  physical  ground  wc  should 
expect  that  the  value  of  m^^(x)  will  be  reduced  by  the  lateral  wake 
because  the  motion  in  the  wake  is  invariably  one  of  low  'total  head' 

(hence  one  of  nearly  uniform  pressure)  and  consequently  its  boundary 
becomes  less  effective  than  a  wetted  hull  surface  in  imparting  kinetic 
energy  to  the  fluid  surrounding  it.  Nevertheless,  theoretical  calculation 


of  the  added-mass  coefficients  of  unsteady  wake  flow  remains  a  most 
difficult  task  as  the  subject  still  needs  major  development.  Of  the  limited 
literature  on  the  subject  reference  can  be  made  to  Wang  &  Wu  (1963)  who 
analyzed  the  small-time  (or  corresponding  to  high-frequency  limit)  wake 
flow  due  to  sudden  acceleration  of  a  flat  plate,  giving  for  the  normal  force 
on  a  broadwise  plate  of  breath  2a  the  result: 

N  =  Nq  +  Nj  +  tN2  +  o(t2)  ,  (105a) 

Nq  =  0.  884  (l+a)(pu2c  a)  ,  (105b) 

Nj  =  m(c)22  Uyc  ,  m(c)22  =  0.  5377  (irPa2)  ,  (105c) 

where  Nq  is  the  normal  force  of  the  basic  flow,  and  a  is  the  wake- 
under-pressure  coefficient.  This  result  shows  that  the  added  mass  of 
the  plate  in  wake  flow  is  0.  53  77  of  its  value  of  a  fully  wetted  plate,  and  is 
very  insensitive  to  a. 

For  an  engineering  solution  to  the  problem  in  hand,  we  recommend 
that  attempt  be  made  to  explore  an  empirical  approach, based  on  the  results 
of  (105  a,  b,  c)  as  a  guide,towards  a  solution  for  moderate  and  low  frequencie 

6.  3  Wave  making  in  shallow  water 

Our  discussion  here  will  be  limited  to  the  range  of  (103)  for  the 
depth  Froude  number,  namely  F^  <  0.  5.  This  range  should  be  sufficient 
for  all  practical  purposes  except  in  unusual  operations.  We  shall  further 
assume  that  the  undisturbed  water  depth,  hQ,  is  uniform. 

One  method  to  calculate  the  wave  field  of  a  ship  moving  in  shallow 
water  is  by  applying  slender-body  theory  and  employing  the  long-wave 
approximation.  Adopting  this  method  Tuck  (1966)  obtained  the  solution 


for  both  subcritical  (F^  <  1)  and  supercritical  (F^  >  1)  cases;  his 
result  for  the  drag  force  for  F^  <  1  is 


rr  s'(x)s'(e)iog  |x-e|dxde ,  ao6) 


2irh  (1-Ft) 
o'  h' 


o  o 


where  S(x)  is  the  cross-sectional  area  of  the  double  model,  S'(x)  =  dS/dx. 
For  the  present  application,  however,  Tuck’s  theory  has  several  drawbacks 
in  that  it  over  simplifies  the  dispersion  effects  in  the  subcritical  range 
the  nonlinear  effects  in  the  transcritical  range  and  in  that  it  will  require 
separate  account  of  the  transverse  current  velocity  and  the  lateral  wake. 

An  alternative  method  is  by  applying  Green's  theorem  as  delineated 
in  Section  6.  This  latter  method  is  expected  to  produce  results  with  higher 
accuracies  over  Tuck's  theory  since  it  only  neglects  the  nonlinear  effects 
which  are  expected  to  be  insignificant  except  for  the  case  of  transcritical 
motions.  It  is  therefore  recommended  that  the  required  Green's  function 
be  first  determined  which  will  also  satisfy  the  kinematic  condition  at  the 
bottom  z  =  -hQ;  a  numerical  scheme  should  be  then  developed,  based  on 
the  analysis. 


8.  Conclusions  and  discussion 

The  present  report  is  Part  II  of  a  series  of  two  reports  wherein  a 
mathematical  model  is  formulated  for  estimating  the  loads  on  a  moored 
vessel  in  an  ocean  current.  In  achieving  this  objective,  it  is  necessary 
to  assemble  a  number  of  different  methodologies  which,  in  total,  constitute 
the  model.  We  believe  that  the  material  developed  herein,  and  in  several 
areas  (such  as  the  treatment  on  the  effect  of  shear  in  separated  flow) 
extends  the  present  state-of-the-art  on  the  subject.  Because  of  this,  it  is 
not  possible  to  forecast  precisely  the  errors  committed  without  actually 
comparing  model  predictions  with  laboratory  and  field  data. 

The  numerous  formulations  and  methodologies  presented  in  the  two 
reports  are  needed  to  cover  the  various  cases,  each  of  which  demands  a 
different  treatment.  At  first,  they  may  appear  somewhat  disjointed.  To 
facilitate  easy  reference,  a  chart  is  prepared  and  shown  in  Table  1  which 
gives  a  roadmap  of  the  overall  structure  and  network  of  the  methodology. 

Referring  to  Table  I,  the  methods  can  be  divided  into  two  main 
groups  in  terms  of  the  angle  of  incidence  of  the  current  with  respect  to 
the  axis  of  the  vessel.  This  is  mandated  by  the  completely  different  fluid 
mechanical  behavior  of  the  flow  field  when  there  is  a  significant  region  of 
separated  flow  such  as  one  that  would  certainly  occur  for  beam-on  current. 
We  believe  the  combination  of  free  streamline  and  strip  theories  recom¬ 
mended  here  for  high  incidence  angles  to  be  the  most  rational  approach. 
While  it  is  not  possible  to  precisely  delineate  its  region  of  applicability 
without  comparing  predictions  with  experimental  observations,  we  believe 
it  would  probably  remain  reasonably  accurate  even  at  incidence  angles 
as  low  as,  say,  30°.  For  small  angles  of  incidence,  low  aspect  ratio 
wing  theory  should  be  used  instead.  We  recommend  computing  the  forces 
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using  both  methods  for  any  angle  less  than,  say  3  0  ,  and  reserving  judge¬ 
ment  as  to  which  one  to  adopt  pending  experimental  verification. 

Among  the  various  formulations  in  the  two  parts  of  the  report, 
those  detailed  in  Part  I  constitute  first  order  models  which  are  expected 
to  provide  quite  good  estimates  of  the  forces  and  moments  in  practice. 

The  formulations  in  Part  II  consists  of  higher  order  models  whose  practical 
effects  would  be  a  modification  to  those  forces  and  moments  predicted  by 
the  methods  in  Part  I.  By  their  nature  as  higher  order  models,  they  tend 
to  be  more  complex  and  thus  more  difficult  to  implement.  Their  contribu¬ 
tions  to  the  augmented  or  improved  forces  and  moments  are  not  expected 
to  be  large. 

It  is  recommended  that  the  methodologies  developed  be  implemented 
numerically  in  a  computational  model  which  can  then  be  used  to  calculate  the 
forces  and  moments  on  a  ship  moored  in  a  current.  Rather  than  constructing 
a  grand  model  incorporating  all  the  formulations  in  these  reports,  it  would 
be  much  more  cost  effective  to  proceed  in  a  stepwise  fashion.  We  recommend 
first  the  implementation  of  the  methods  in  Part  I  in  a  computational  model. 

It  should  then  be  used  to  predict  the  forces  and  moments  on  a  number  of 
generic  ship  hullshapes  such  as  destroyers,  tankers,  and  cruisers.  The 
predictions  can  then  be  examined  in  light  of  experimental  findings,  both 
from  the  laboratory  and  the  field.  A  rational  decision  can  then  be  reached 
as  to  the  choice  of  further  implementation. 

The  ultimate  practical  engineering  goal  is  the  estimation  of  current 
loads  on  moored  vessels.  It  would  be  highly  desirable  to  be  able  to  obtain 
rapid  and  reasonably  accurate  estimates  on  the  basis  of  a  relatively  small 
number  of  parameters  characterizing  the  ship  (such  as  length,  beam,  draft, 
block  coefficient,  ship  type,  etc.  )  rather  than  detailed  shiphull  cross 
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sections  along  the  full  length  of  the  ship.  We  foresee  the  possibility  that 
results  from  calculations  using  the  detailed  model  and  generic  ship  hull- 
shapes,  when  properly  interpreted  in  conjunction  with  laboratory  and  field 
observations,  might  be  utilized  to  form  the  basis  of  an  engineering  model. 

A  justification  of  the  simplifying  assumption  of  neglecting  the 

viscous  effects  in  treating  the  wake  flow  in  a  sheared  stream  can  be 

provided  as  follows.  There  are  two  roles  the  viscous  effects  can  play  in 

this  case.  One  is  giving  rise  to  a  boundary  layer  over  the  wetted  part  of 

the  hull  (upstream  of  the  wake  region),  which  is  thin  compared  to  the  ship 

dimension.  The  other  is  to  increase  transverse  diffusion  of  the  vorticity 

1  /2 

at  the  rate  of  (vt)  ,  where  v  is  the  kinematic  viscosity  of  water  and  t 
is  the  time.  This  length  of  diffusion  would  amount  to  less  than  1  cm  in  the 
time  when  a  current  of  1  knot  has  traversed  the  ship's  beam  of  15m.  The 
impact  of  viscous  diffusion  on  altering  the  velocity  of  a  rotational  flow  is 
therefore  insignificant  and  can  be  neglected. 

A  final  comment  is  in  order  concerning  the  applicability  of  the 
present  model  for  the  case  of  a  random  wave  field.  Throughout  the 
development  of  both  the  original  model  and  the  present  extensions,  no 
incident  wave  field  is  considered  since  it  is  not  part  of  the  present  contract 
work  statement.  The  wave  radiation  resistance  consideration  described  in 
section  6  pertains  only  to  the  wave  making  as  a  result  of  the  interaction  of 
the  current  with  the  ship  and  cannot  be  extened  to  apply  to  the  effect  of  an 
ambient  wave  field,  whether  it  is  regular  or  random.  The  latter  problem 
is  one  of  wave  scattering.  Techniques  are  available  for  calculating  the 
flow  field  due  to  incident  waves  scattered  and  diffracted  by  a  ship  either 
stationary  or  self-propelled  in  the  forward  direction,  but  in  either  case 
without  a  transverse  current.  In  such  cases  the  incident  waves  are  usually 
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taken  to  be  monochromatic  and  applications  to  the  case  of  a  random  sea 
are  made  by  means  of  linear  superposition  based  on  a  known  wave  energy 
spectrum.  However,  in  asmuch  as  the  present  problem  is  involved  with 
the  possible  presence  of  a  transverse  wake,  the  wave  scattering  problem 
would  require  new  techniques  for  its  solution.  At  present,  there  is  a 
virtually  complete  lack  of  information  for  this  case,  in  either  theoretical 
analyses,  or  experimental  observations. 

Since  the  present  problem  is  aimed  at  a  moored  ship  and  hence  in 
a  rather  sheltered  water,  incident  waves,  if  any,  are  not  expected  to  be 
large  and  their  forces  would  be  small  compared  to  the  other  hydrodynamic 
forces  considered  herein.  In  the  unexpected  circumstance  of  a  sudden 
storm  resulting  in  large  incident  waves,  such  as  when  the  peak  orbital 
velocity  of  the  water  particles  is  equal  to  the  mean  current  speed,  the 
combined  forces  will  be  so  strongly  time  dependent  that  the  quasi- steady 
theory  will  no  longer  be  adequate. 
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